Photoassociation dynamics in a Bose-Einstein condensate 
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A dynamical many body theory of single color photoassociation in a Bose-Einstein condensate is presented. 
The theory describes the time evolution of a condensed atomic ensemble under the influence of an arbitrarily 
varying near resonant laser pulse, which strongly modifies the binary scattering properties. In particular, when 
considering situations with rapid variations and high light intensities the approach described in this article leads, 
in a consistent way, beyond standard mean field techniques. This allows to address the question of limits to the 
photoassociation rate due to many body effects which has caused extensive discussions in the recent past. Both, 
the possible loss rate of condensate atoms and the amount of stable ground state molecules achievable within a 
certain time are found to be stronger limited than according to mean field theory. By systematically treating the 
dynamics of the connected Green's function for pair correlations the resonantly driven population of the excited 
molecular state as well as scattering into the continuum of non-condensed atomic states are taken into account. 
A detailed analysis of the low energy stationary scattering properties of two atoms modified by the near resonant 
photoassociation laser, in particular of the dressed state spectrum of the relative motion prepares for the analysis 
of the many body dynamics. The consequences of the finite lifetime of the resonantly coupled bound state are 
discussed in the two body as well as in the many body context. Extending the two body description to scattering 
in a tight trap reveals the modifications to the near resonant adiabatic dressed levels caused by the decay of the 
excited molecular state. 
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I. INTRODUCTION 

Since the pioneering achievement of Bose-Einstein conden- 
sation in dilute alkali gases much effort has been made to ex- 
tend this success to a wider class of atomic and molecular 
species. In particular, the production of quantum degener- 
ate ensembles of molecules 1 1 ] promises a wide spectrum of 
important applications ranging from ultra-precise molecular 
spectroscopy and cold collision studies Ql to "superchemical" 
reactions 1 3 ] and the investigation of a possible crossover from 
a Bose-Einstein condensate (BEC) of molecules to Bardeen- 
Cooper-Schrieffer (BCS) correlated pairs in Fermi gases 0]. 
Laser cooling of molecules is essentially ruled out due to the 
complicated rovibrational molecular level structure | y], such 
that alternative methods are in order when proceeding towards 
degenerate molecular gases. The most promising procedure 
to date is photoassociation (PA) in ultracold, possibly degen- 
erate atomic gases |2, 5]. A photoassociative scattering res- 
onance may occur when a conveniently tuned laser field in- 
duces a transition of a colliding atom pair to an electronically 
excited molecular, i.e. an excimer bound state. PA in BECs 
has first been reported in Ref. |6, 7] and recently been per- 
formed at considerably higher rates 1 8, 90. Here, we present a 
dynamical many body theory which generalizes the approach 
in Refs. | \^ [Tl], [l31 to photoassociation in a BEC and sys- 
tematically takes into account effects of correlation functions 
beyond the condensate mean field. 

On the broadway towards molecular BEC considerable suc- 
cess has been achieved recently in a number of exciting exper- 
iments employing Feshbach resonances. Such resonances are 
formally analogous to PA resonances but they are not compli- 



*email:T. Gasenzer@thphys.uni-heidelberg.de 



cated by a short lifetime of the resonantly coupled molecular 
state. Macroscopic amounts of ultracold diatomic molecules 
were assembled both from Bose-Einstein condensed atoms 
llofl4iri5iriq. flrl and in degenerate two-component gases of 
fermionic atoms 1 18, 19, 20, 21, 22, 23, 24, 25. 26, 27, 28, 29]. 

Despite these remarkable achievements PA in the degener- 
ate regime remains at the top of the list of desirable aims. PA 
promises that one may use far developed laser techniques to 
form molecular quantum gases of arbitrarily modelled com- 
pounds. Techniques available so far are limited due to insuffi- 
cient overlap between molecular states in different scattering 
channels as well as by the rapid spontaneous decay of the in- 
termediate excimer states. This decay may be controlled by 
applying a second laser beam which drives the transition of the 
excited molecules to a stable bound state of atoms in the elec- 
tronic ground state |3, 5, 6], particularly in the form of Stim- 
ulated Raman Adiabatic Passage (STIRAP) 03 HH. Coher- 
ent control techniques which achieve optimized laser pulses 
II32I1 may further help to produce condensates of molecules 
in conveniently engineered states. Theoretically, PA involv- 
ing one or more laser modes has been studied in detail in 
the framework of two body scattering 1 2lll l33ll34ll3ll3^ 
as well as for quantum degenerate gases containing bosonic 
L37..38..39..40..41.1 and/or fermionic atoms |42|. An important 
question which remains to be answered both experimentally 
and theoretically concerns limitations to achievable associa- 
tion rates due to many body effects 18, 9, 36, 39, 40, 41. 431. 

Early approaches to the many body dynamics of a BEC near 
a photoassociation resonance were based on mean field theory 
in the form of coupled Gross-Pitaevskii equations for atoms 
and molecules ll44ll . Beyond mean field level the effect of the 
pai r correlation function has been taken into account 1 11, w, 
|4fl , 47]. All of these approaches, including pure mean field 
models, are expected to satifactorily describe condensate loss 
and molecule formation in adiabatic crossings of a Feshbach 



resonance, in which the magnetic field is varied linearly in 
time 1 45, 48]. It was shown, however, in Ref. LI 1.1 . that a non- 
trivial variation of the magnetic field near such a resonance 
may lead to effects which are not accounted for by coupled 
Gross-Pitaevskii equations. 

In view of these results a systematic treatment of the dy- 
namics of higher order correlation functions seems is in order, 
which allows to describe strong coupling PA in a BEC and to 
identify phenomena which are distinctly different from those 
near a Feshbach resonance. 

The article is structured as follows: Section ITTI is focused 
on the two body physics of low energetic atoms colliding in 
the presence of the near resonant laser mode. The dynami- 
cal evolution of the many body system starting as a trapped 
Bose-Einstein condensate is studied in Section|lIl] Section lTVI 
contains a summary of the results and draws our conclusions. 
An appendix contains explicit descriptions of the techniques 
used to describe the two body collisions and the many body 
dynamics as relevant for this work. 

In Section |II] we define the two-channel picture of reso- 
nant single color photoassociative scattering of a pair of ^'^Na 
atoms as in the recent experiment of McKenzie et al. ^. The 
scattering length approach to the laser assisted atomic colli- 
sions |33, 3?] is introduced, and the relevant parameters char- 
acterizing the resonance are identified. These parameters de- 
fine a minimal model Hamiltonian which suffices to describe 
all the relevant low energy scattering properties. Using this 
Hamiltonian we focus on the set of atomic pair states dressed 
by the laser with a contribution from the excimer bound state. 
To illustrate the effects of the finite lifetime of this bound state 
the dressed state spectrum is studied for pairs confined in a 
tight harmonic trap for which a convenient modification of 
the minimal model Hamiltonian is developed. In this way we 
could study the existence of quasi-bound states above thresh- 
old which are a consequence of the excimer state decay. We 
then consider the components of the dressed state closest to 
threshold corresponding to the background channel and the 
admixed excimer state. We show that for high light intensities 
the near resonant dressed state resembles a bound state in the 
background channel which, as discussed in Section IHl] has 
a decisive impact on the corresponding many body dynamics 
and on possible single color PA rates. 

In Section Hn] we extend the many body dynamical theory 
developed in Refs. Olllllllillll to single color PA in 
a BEC. We discuss in detail the incorporation of the resonant 
two body physics into the many body formalism and the role 
of effects beyond the Gross-Pitaevskii mean field level. Since 
the theory exactly respects atom number conservation at all 
times it also allows us to study the overall loss of atoms due 
to the spontaneous decay of the excited molecular state. As a 
result we obtain limits to both, the achievable condensate loss 
rate and the maximum rate at which ground state molecules 
may be formed by spontaneous electronic deexcitation. The 
different origins of these limits will be discussed. Our re- 
sults show that the rate limit for condensate loss predicted in 
Ref. 1 41] is exceeded at high light intensities. Nevertheless, 
we find that the maximum amount of atoms lost within a cer- 
tain time interval is reduced compared to the value derived 
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FIG. 1: Schematic interaction potentials of ^^Na pairs coupled by a 
single laser mode. The condensate atoms collide with a relative en- 
ergy close to the background channel threshold indicated by a dashed 
line. They are excited by a laser with frequency u) — {2n) 80.7 THz 
to the (J = Ijf = 135) rovibrational bound state in the closed 
channel and spontaneously decay by photon emission to ground 
state molecules as well as free atom pairs, with a time constant of 
7"^ = 8.6 ns. A = 



KeE^/h denotes the detuning of the laser. 



from a naive Gross-Pitaevskii approach. This limit, however 
has not yet been reached at the intensities in the experiment in 
Ref. ]8]. We finally present examples for the time evolution 
of the density profile of a trapped ensemble under high light 
intensity PA which exhibit strong collective excitations of the 
atomic cloud. 



II. LASER ASSISTED RESONANT ATOMIC COLLISIONS 

Consider two atoms which are colliding in the presence of 
a light field from a laser source. PA of the pair is generally un- 
derstood as its light induced binding to a diatomic molecule. 
We would like to choose a somewhat more general definition 
which includes any laser controlled modification of the colli- 
sional interaction properties of the atoms. A key finding then 
is that PA allows to modify the low energy scattering proper- 
ties of atoms in an ultracold vapour in a way very similar to 
what is known as a Feshbach resonance I33ll34ll . Although 
very different physically much of the theoretical formalism 
used for one kind of scattering resonance, photoassociative or 
Feshbach resonant, may also be applied to the other. 

The main aim of the work presented in this article is to adapt 
the microscopic dynamical many body theory of Refs. llOlfllll 
to photoassociation of Bose-Einstein condensed atoms. In this 
theory the binary interactions are incorporated in a form in 
which they are treated in quantum mechanical scattering the- 
ory. We shall therefore introduce, in this section, the quantum 
mechanical description of the low energy binary collisions in 
the presence of the laser field. The photoassociative colli- 
sions are modelled in a two channel picture which accounts 
for asyptotically free atom pairs with a certain relative en- 
ergy coupled by the laser to a closed channel molecular bound 
state, cf. Fig.^ We follow the approach of Refs. L33..34.1 in 
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FIG. 2: Complex s-wave scattering length a(/, A) describing the 
low energy scattering of "^^Na atoms, in the vicinity of a scatter- 
ing resonance induced by single color photoassociation, for / — 
1 kW/cm^, as a function of the detuning A (see Fig.0. 



characterizing the low energy scattering in terms of a laser 
modified s-wave scattering length. Due to its relatively strong 
dipole moment the closed channel state is short lived such that 
the physical situation is substantially modified as compared to 
the Feshbach case. We shall study in detail the consequences 
of this decay for the spectrum of dressed states of an atom pair 
in the vicinity of the resonance. 



A. Photoassociation of Na atoms 

In this article we focus on single color photoassociation in 
a trapped condensate of ^^Na atoms. Related experimental 
results have been presented by McKenzie et al. in Ref. |^. 
In this experiment the condensed atoms are initially in the 
3^Si/2, \F = l,?7ii? = —1) hyperfine state, and their bi- 
nary collisions in the absence of laser light are described by 
a mixture of the X^E+ singlet and the a'^I]+ triplet potentials 
|3l- Atom pairs in the condensate which initially are mov- 
ing freely on the nearly degenerate asymptote of these ground 
state potentials are coupled by a laser with frequency oj — 
(27r)80.7THz to the (J = 1, i^ = 135) rovibrational level 
of the A^E+ Na2 molecular potential (Fig.^. This excited 
molecular state |0i,) has a lifetime of about 8.6 ns which cor- 
responds to a linewidth of 7 = (27r) 18.5 MHz being nearly 
twice the atomic linewidth. The spontaneous decay of |0y) to 
hot atom pairs or ground state molecules is taken into account 
in the imaginary part of the energy E^ ~ hiuj^ — 17/2) of 
\(\>v)- If the intensity of the laser light is kept fixed, the tun- 
ing of the frequency through the transition leads to a resonant 
variation of the s-wave scattering length 133. 34.i36ll 



a(/, A) = flbg 1 



m 



(1) 



Ao(/) - A - ^7/2 

as shown in Fig. |2] The description of low energy PA scat- 
tering in terms of a scattering length indicates the formal sim- 
ilarity to a Feshbach resonance. However, due to the finite 



lifetime of |0y) the PA scattering length is complex and its 
modulus does not diverge on resonance. A = w — o^i, is the 
detuning of the laser from the frequency difference between 
the molecular level v and the ground state asymptote. The 
resonance is shifted by Ao(/) — ojo{I) — uj^. On resonance, 
uj — uJo{I), the real part of the scattering length is equal to 
the background scattering length abg which characterizes the 
interactions in absence of the laser. The width of the reso- 
nance is parametrized by S{I). Both, (5(/) and the shift Ao(/) 
depend on the light intensity / and, for low intensities, were 
found to vary linearly with / @] as predicted by perturbation 
theory |36]. The scattering length shown in Fig.|2]was deter- 
mined with Eq. Q using the values 6{I) = (27r) 38.4 MHz, 
and Ao(/) = (27r) 164 MHz, for I ^ l.OkW/cm^ 0]. 



B. Two channel model of resonant two body scattering 

Figure ^ shows the general process we aim to describe. 
Two ultracold condensate atoms with a relative kinetic energy 
close to the free-bound threshold of the potential Vbg interact 
with each other In the absence of the laser light their colli- 
sional interaction is described by Vbg- Since the total energy 
of the atom pair lies above threshold this background scat- 
tering channel is termed an open channel. For simplicity we 
consider only a single open channel potential, e.g. , in the case 
of sodium atoms colliding in the \F = l,mF = — 1) state, a 
potential which reproduces the low energy scattering proper- 
ties of condensed sodium atoms as measured in experiment 
or calculated theoretically using more elaborate multichannel 
techniques |3l • We neglect the explicit notation of the internal 
state, of fine and hyperfine quantum numbers, in order to fo- 
cus on the structure of the light induced coupling of the states. 

The potential which describes the interactions of atoms in 
the Si/2 ground state with atoms in the excited P1/2 state is 
denoted as Vd- The difference of the threshold energies of Vc\ 
and Vbg equals the atomic line frequency and, since this is in 
practice larger than uj, Vd is termed a closed channel potential 
for the near zero energy condensate atoms. 

We assume the laser to be described as a classical plane 
wave with wave vector k and frequency uj = |k|c. The 
light gives rise to a radiative coupling W(x, i) = hil{x,t) 
expressed in terms of a time dependent Rabi frequency 
ri(x, t) ~ rJo(t) cos[kx ~ ujt]. Here, we have included the 
temporal variation of the envelope of the light field, i.e. the 
time dependence of its intensity /(<) in the Rabi frequency 

noit). 

The Hamiltonian operator describing the binary interac- 
tions in the two channel picture may now be written as H2B = 
Ho + V2B, with 



Ho^ rkin[|bg)(bg| 

V2B= (|bg),|cl)) 



|cl)(cl| 

Vbg w 
W Vci 



(bg 
(cl 



(2) 
(3) 



Here, Tkin = —fi^W'^/m is the relative kinetic energy. |bg) 
and |cl) are the superposition of products of atomic states as- 
sociated with the open and closed channels, respectively. In 



Eq. (|2ji the trapping potential may be neglected since it does 
not vary on the scale of the binary interactions. 

The range of relative kinetic energies in ultracold collisions 
is small compared to the typical spacing between vibrational 
states. Throughout this work we will therefore consider only 
situations where the laser couples to a single well defined 
closed channel bound state \(t)u)- The closed channel Hamil- 
tonian H^x — Tkin + Vc\ can then be replaced by the operator 
E^Py, with the quasi projector P^, = \4>v){4'u\ onto the ex- 
cimer state, Hc\\(t>v) = Ey\(j)y). Because of the spontaneous 
decay of the 3^Pi/2 excited atoms the ((/)y | are left eigenstates 
of the non-Hermitian operator Hc\ to the same complex eigen- 
value Ey, {(j)u\Hc\ = {(t),y\Ey. Cf. AppendixiBlfor details. 

Within the subspace spanned by the open channel states and 
the single closed channel state 10^) the Hamiltonian reads: 

Ho= rkin|bg)(bg| + |0,,cl).B,(0,,cl| (4) 

(b 



v.B^{H)M))[^^^'^^' 



(5) 



For the following discussion of the dressed state spectrum of 
the Hamiltonian it will be convenient to go to an interaction 
picture, where, in the rotating wave approximation, the phase 
exp[—iujt] of W{x, t) is removed and the closed channel fre- 
quency w^ = Reii'^/^iinEq. (0}isreplacedby — A = oj^ — oj. 



C. Two channel dressed states 

It is well known that the energy of the uppermost vibra- 
tional bound state of two interacting particles determines the 
singularities of the s-wave scattering length. 

A Feshbach resonance may be induced through Zeeman 
shifting atomic hyperfine levels and thus different scattering 
channels relative to each other by applying an external mag- 
netic field. The channel states are coupled to each other 
through internal atomic forces. This leads to a modification 
of the bound state spectrum of the atom pair In contrast to 
this, the external electromagnetic field of a PA laser transfers 
energy to the atom pair. As a result, the width and the shift of 
the scattering resonance depend on the intensity of the light. 

The two body state resulting from the modification by the 
laser field is called a dressed state. Moreover, since the life- 
time of the closed channel state is very short the population 
of the dressed states may not be measured in the same way 
as that of the bound states in recent Feshbach crossing exper- 
iments L15..16..18..19..20..2L.25.1 . viz. by reversing the asso- 
ciation process and imaging the atoms. Stable molecules are 
rather obtained as bound states of the ground state potential 
Vbg when the closed channel molecule spontaneously emits a 
photon or when the transition is induced, e.g., by illuminating 
the sample with one or more additional lasers. Such processes 
of Raman PA 1 5 ] , possibly in the form of STIRAP 1 3 1 1 or with 
optimized laser pulses (Coherent Control) 113211 are beyond the 
scope of this article, but may be described in a straightforward 
extension of the formalism presented here. 

Nevertheless, the dressed state picture is formally equiva- 
lent to the description of the modified bound state spectrum 



near a Feshbach resonance. In the interaction picture the 'en- 
ergy' of the closed channel state compared to the open chan- 
nel threshold is given by the laser detuning which one may, 
e.g., ramp through the resonance in much the same way as 
it is done in the case of adiabatic molecule association at a 
Feshbach resonance. 

For these reasons we shall in this subsection discuss in de- 
tail the spectrum of dressed states of interacting atoms pairs 
in the framework of the two-channel picture introduced in 
Sect. Ill B1 After a discussion of the low energy scattering in 
the background channel we will derive the continuum of scat- 
tering states as well as the discrete part of the dressed state 
spectrum close to the PA resonance. In Section lTlDl we will 
describe the particular consequences of the spontaneous de- 
cay of the closed channel state for the structure of the dressed 
state spectrum. 



I. Background scattering 

As was discussed in detail in Refs. Illl |451 the binding 
energy of the highest vibrational state of Vbg is determined, 
to an excellent approximation, by Obg and the long range 
r-dependence of Vhg{r). For the case of ground state al- 
kali atoms, this is dominated by the van der Waals inter- 
action, Vbg(r) ~ — Cg/r^, setting a length scale, Zvdw = 
(mCg/Sfi.^)^/^ which is much smaller than the thermal de 
Broglie wave length in typical ultracold collisions. Using 
semiclassical arguments it was shown in Ref. 15511 that the 
first correction to the energy of the highest vibrational state 
is given through 



E-i 



m{ahg - ay 



(6) 



i.e. by the usual formula for the energy in terms of the 
scattering length which, however, is shifted due to the 
van der Waals tail by the mean scattering length a — 
Uwr(3/4)/(V2r(5/4)) ~ 0.956 /vdw- For the case of 
sodium considered here one has Cg — 1561 a.u. ll64ll (1 a.u. = 

0.095734-10"^'^ J nm'5) and thus /vdw = 45.0aBohr (aBohr = 
0.052918 nm), and, with Obg = 54.6 OBohr pS], a binding en- 
ergy of S_i = (27r) 1164 MHz ft. 

The scattering of sodium atoms in the |F = 1, nip = — 1) 
state is described by a mixture of the X^I]+ singlet and a''l]+ 
triplet potentials with scattering lengths for which Ref. I58ll 
gives values of a^^o = 20.3aBohr and as=i = 63.9aBohr, 
respectively. The energies of the last bound states of the 
singlet and triplet potentials have been experimentally deter- 
mined as E-i{S = 0) = -318(21) MHz ft and E-i{S = 
1) = -270(36) MHz ft, respectively |59. 60]. Hence the for- 
mula ^ is not applicable for the sodium potentials. This is 
due to the fact that the binding energies are too large such that 
the mean scattering length a = 43 OBohr is not a small cor- 
rection to the background scattering length abg = 54.6 asohr 
155,56]. 

However, it was shown in Ref. |57] that the energies of the 
uppermost bound state of the singlet and triplet potentials are 



excellently determined by the scattering length and the van 
der Waals coefficient Cg alone. At zero collision energy the 
scattering length Obg summarizes all the unresolved details of 
Vhg{r) into a single length scale. Atoms interacting at finite 
collision energies resolve details of the potential, and the long 
range behaviour quantified by the van der Waals length deter- 
mines the leading contribution. Any model of the background 
scattering potential recovers the binding energy of the highest 
vibrational state and the low energy scattering properties if it 
properly accounts for the parameters Cbg and Cg . 

In this work we shall use the experimentally determined 
value of E^i = — 318MHz h to model the long range be- 
haviour of the background potential. We introduce a minimal 
background potential in Appendix|X]which is related to a sep- 
arable approximation of the fuU interaction potential. 



2. Dressed state spectrum: Continuum part 

Assuming that the spectrum and states associated with the 
background potential are known, we shall, in this section, de- 
rive the continuum of scattering states of the full two chan- 
nel Hamiltonian, Eqs. (|4}, (|5}- Moreover we will deduce the 
dependence of the scattering length a{I, cu) on the relevant 
quantities contained in the Hamiltonian. This will enable us 
to describe the full resonant two channel scattering properties 
at low collision energies by a minimal set of five parameters. 

The two channel eigenstates of the Hamiltonian, Eqs. 0, 
(|5}, may be written in terms of their open and closed channel 
components as |bg) jc/)'^*;^ ^ |cl) |(/)'^'). The two components are 
solutions of the stationary two channel Schrodinger equation, 
which, in the interaction picture introduced in Sect. Ill Bl reads 






w 



= E 



(7) 



where iJbg = Tkjn + Vbg andi^ci = {E^~hLj)P^ = -h{A + 
i7/2)P„withP„ = |0,)(0,|. 

As discussed in Sect. Ill B I the laser is assumed to couple 
the ground state atom pair to a single unstable bound state in 
the closed channel potential. As a consequence, the closed 
channel component of the scattering state vanishes at large 
interatomic distances. 



6p(r) 



0. 



(8) 



Assuming that there is no other decay channel, the continuum 
of scattering states then consists of all solutions of Eq. with 
real positive energy E = p^ /m, corresponding to a relative 
momentum p ~ |p| of the atoms colliding in the open chan- 
nel. An imaginary part of E is excluded since the atoms are 
in the ground state and thus stable at large distances r where 
the coupling to the closed channel resonance state is negligi- 
ble. The background scattering states at large r are chosen to 
obey the well known expression in terms of an incoming plane 
wave, an outgoing spherical wave, and a scattering amplitude 
/(i9,p) which depends on p and the angle d between p and 



the relative position r of the outgoing atoms. 



cpl^iv) ^ i2nh)-^/' 



ipr/h 



+ f{^,p)- 



(9) 



The Schrodinger equation Q may then be expressed in 
terms of the Green's functions 



Ghgiz) = 
Gei(z,A) = 



(^-^bg) \ 

Pu 

z + h(A + 17/2) ' 



(10) 
(11) 



which depend on a complex energy parameter z 14911 . For the 
scattering states the equation reads: 






Gbg(^ + iO)W^P,|0^') 
Gci{E,A)P,W\4>'^s) 



(12) 



The imaginary part 'iO' of the argument of Gbg indicates that 
the Green's function is evaluated in the limit of approaching 
the energy E from the upper half of the complex plane, ensur- 
ing that the wave function (j>p^{r) has the asymptotic behavior 
given in Eq. (|5). 

In Eq. (I12> . |0p ) denotes the positive energy solution of 
the homogeneous Schrodinger equation, i.e. Eq. with W 
set to zero, behaving at large interatomic distances as 



0W(r) 



(27rft)- 



-3/2 



Jpr/h 



./bg(^,rt- 



oV" 



(13) 



Given the pole approximated closed channel Green's func- 
tion Gc\{E, A), Eq. din , the closed channel component of the 
scattering state, Eq. (I12t . is proportional to the resonant state 
101/). Thus the scattering state may be written as 



l/p' 







(14) 



with the amplitude 

AiE,A) 



K\W\ 



j^bgN 



E + h{A + ij/2) 



(15) 



which from the background channel component of Eq. ( I14t . 
is determined as 



AiE,A)^ 



{4>.\W\^i+^] 



E + h{A + t-f/2) - {cP^\WGbg{E + iQ)W\<j>,,) ' 

(16) 



In summary: Given the full positive energy spectrum 1 0p ) 
of eigenstates of i/bg as well as the resonantly coupled closed 
channel state \(j),^), Eqs. ( I14> and (I16> represent the exact so- 
lution of the two channel Schrodinger equation 0. Under the 
assumption that the excitation of the colliding atom pair is re- 
stricted to the state \(f),y) this pole approximated solution of the 
full two channel Schrodinger equation becomes exact. 



We shall now derive formula Q for the scattering length 
near resonance. The s-wave scattering length in the open scat- 
tering channel is defined as the negative of the zero momen- 
tum limit of the scattering amplitude f{i3,p) which is equiva- 
lent to the expression 

a =. lun{2nh)^-—{p,hg\T2B{pVm)\p,hg) (17) 



p^O 



Airff 



in terms of the on shell T-matrix element between plane wave 
states with momentum p in the background channel. Using 
Eqs. 0, (inj, (|T6}, and the definition (p|r2B(p^/m)|p) = 
(p|V2b|0p'^'^) of the diagonal plane wave T-matrix elements 
in terms of the potential V2B, Eq. (|5}, and a two channel scat- 
tering state l^p'^'') as given in Eq. (I14> . we find: 



{2TrhY 



ibg 



«bs 






(18) 



This reproduces the resonance formula for the scattering 
length given in Eq. Q if one identifies 



S{I) 



mh 



{2nh)^ct>[i+^\W{I)P.W{I)\4+^}, (19) 

(20) 



47r?i2abg 
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where we explicitly denote the dependence on the light inten- 
sity /. 

Knowledge of the dipole matrix elements and the Franck- 
Condon factors would allow, on the basis of the known back- 
ground scattering properties, to calculate the shift Aq and 
widths of the photoassociation resonance from the matrix el- 
ements in Eqs. ( I19> , ( I20> (cf., e.g. Ref. I6lll \ In this work 
we will however assume that the shift and width are quantities 
known from experiment over a certain range of light intensi- 
ties /, and that they depend, in this range, linearly on /. The 
experiment in Ref. 1 8] confirms this linear dependence for the 
range of applicable intensities / applied, and we will assume 
that the linear behaviour also extends to somewhat larger in- 
tensities. 

In analogy to the minimal model potential describing the 
low energy scattering properties of the background channel 
interactions we shall introduce a characterization of the cou- 
pling to the closed channel state in terms of two parameters. 
These two parameters measure the coupling strength and the 
range of interatomic distances where the coupling occurs and 
are fixed through Eqs. ( I19> and ( I20> in terms of the measure- 
able resonance width and shift. Cf. Appendix IaI for details 
of how to choose an ansatz convenient for practical computa- 
tions. 

In summary, together with the two background parameters 
of the minimal model, and the decay width 7 of \(f>i,), they 
form a minimal set of five parameters which fully describe 
the low energy two body scattering of atoms under the influ- 
ence of the near resonant photoassociation laser In giving 
explicit expressions for matrix elements of the Hamiltonian 
operator this minimal parametrization represents a powerful 
tool when solving the many body dynamical equations which 
will be presented in Sectionlim 



3. Dressed state spectrum: Discrete part 

The discrete part of the spectrum of the two body Hamil- 
tonian consists of dressed states with in general complex 
eigenenergies since the corresponding spatial wave functions 
vanish at large interatomic distances. Similarly to Eq. (I12> the 
Schrodinger equation may be written in integral form. 
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and has, taking into account Eq. (Ill> . the normalized solution: 
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\(t>J) has by assumption unit norm such that the overall nor- 
malization constant is 

ACi - [1 + (0.|M/^Gbg(£;d)^Gbg(£;d)W^|0.)]'^' , (23) 

and the complex dressed state energy is determined by the 
condition 

-Ed = {4>AWG^^{Ea)W\(I>,) - h{A + Z7/2). (24) 

The states given through Eqs. ( I22t and ( I24> are right eigen- 
states of the non-Hermitian Hamiltonian H2B, which are, in 
general, not pairwise orthogonal to each other and to the scat- 
tering states given in Eq. (I14> (Cf. AppendixiBl. 



D. Properties of the dressed states near resonance 

1. Adiabatic energy levels: homogeneous case 

Suppose, for a moment, that 7 = 0. In this case Eq. MA\ 
recovers relation ( I20t when setting E^ — Q and taking into 
account that A is tuned to resonance, A = Ao(/). This for- 
mally corresponds to the situation at a Feshbach resonance 
and illustrates that the scattering length diverges where the 
dressed state energy vanishes. The adiabatic energy level 
structure close to the scattering resonance, where the closed 
state crosses the background channel threshold, is shown in 
Fig. The graph shows the dependence of the then real 
dressed state energy E^ on the detuning A of the laser from 
the frequency of the transition. 

We calculated this dependence using Eq. ( I24> for photoas- 
sociation of ^■^Na atoms described in Sect. Ill Al but with 7 set 
to zero. Technically this is achieved using the minimal model 
Hamiltonian which is described in Appendix |X] and fixed by 
the background scattering length Obg = 54.9 aeohr (aBohr = 
0.052918nm) 158], the binding energy £;_i//i = -318MHz 
JH, the width d6{l)/dl = 38.4MHz/i/(kW cm'^), the 
shift dAoiI)/dI = -164MHz/i/(kW cm'^) 0], and the 
decay width 7 = 0. 

When approaching the resonance from the side where the 
scattering length is positive the bound state energy E^ is 
negative and, close to resonance, universally behaves as in 
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FIG. 3: Energy Ed of the uppermost bound state of two sodium 
atoms for the hypothetical case of a stable closed channel bound 
state, 7 = (thick solid line, for A < Ao: — ■ — line). Ed 
is shown as a function of the detuning Ao(/) — A from the reso- 
nance, for a light intensity I — IkWcm"^. For large detunings 
the bound state is similar to the last bound state of the background 
potential (thin solid line). The inset shows _Bd(Ao — A) close to 

resonance, where it is approximated by Ed = —h?/{rna?) ( — 

line) or Ed = —t^/{m[a — a^) (dashed line), see Sect. lIICTI Ap- 
proaching resonance (A -^ Ao(/)), the scattering length diverges 
and Ed -^ 0. 
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Eq. (|6}, with a hp^ r eplaced by the near resonant scattering 
length, Eq. Q |451. This dependence of the dressed state 
energy on A is modified when the closed channel state has 
a nonzero width 7 = (27r) 18.5 MHz. An example of such 
a modified behaviour is shown in Fig.|3 for a light intensity 
/ = lOkW/cm^. In general, E^ is complex, and at the reso- 
nance, A = Ao (/), neither its real nor its imaginary part van- 
ish. The real part rather crosses the threshold with a non-zero 
slope and becomes positive at Aq — A ~ — (27r) 18.73 MHz. 
At Aq — A ~ (27r) 0.12 MHz the imaginary part vanishes and 
the dressed state becomes degenerate with a scattering state 
with non-zero energy. 



We may estimate this energy by considering the analytic 
behaviour of E^ near resonance. As in the case of a Feshbach 
resonance, where the binding energy approaches E\^{B) = 
— h? / m[a{B)]'^ when a{B) -^ +00 [45], the same uni- 
versal dependence on a applies to PA if \a\ is sufficiently 
large. In Fig. 0] we compare the real and imaginary parts 
of the energy given by Ed{I,A) = — fi,^/m[a(/, A)]^ and 
Ed{I,A) = -h'^/m[a{I,A) - a]^ where a{I,A) is deter- 
mined from Eq. Q, with the energy derived from Eq. ( I24> . 
When crossing the resonance from the negative detuning side 
the dressed state is therefore transferred smoothly into a scat- 
tering state with energy E ~ h'^'^^/{Ama^ S^). This esti- 
mate is valid if S ^ 7/4. For / — lOkW/cm^ one finds 



E =ReEd,ma.x/h = 0.12 MHz. We shall study this intrusion 
of the dressed state into the positive imaginary plane in more 
detail in the following subsection. 



FIG. 4: Energies Ed of the discrete dressed states closest to thresh- 
old, as a function of Ao(/) — A. The graph shows the same levels as 
Fig-El but for a light intensity of / = 10 kWcm^^ and a decay width 
7 — {2tt) 18.5 MHz which renders the level energies complex. The 
dotted and dashed lines show the analytic approximations to the near 
resonant behaviour as in Fig. |3| Note that the real part of -Ed may 
become positive, i.e. intrude the continuum of scattering levels with 
real energies Ep > 0. 



2. Adiabatic energy levels in a tight trap 

In Sect. Ill DTI we have shown that the real part of the energy 
of the near threshold discrete dressed state may become posi- 
tive for laser detunings near resonance. Furthermore, we have 
seen, that at the largest possible real part the imaginary part 
vanishes such that the state becomes degenerate with a state 
of the continuum of stable scattering states. In the following 
we shall illustrate this behaviour by considering atom pairs 
that are confined in their motion by a tight spherical harmonic 
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For the first we consider again a stable molecular state, 
7 = 0, such that the Hamiltonian is hermitian. The adia- 
batic levels, i.e. the dependence of the energies on the detun- 
ing A is generally characterized by avoided crossings of the 
diabatic open and closed channel states as shown in Fig. |5] 
(cf. also, e.g., Ref. I6al '). The energies in this graph were 
computed for the same parameters as in Fig. |3j with the ad- 
dition of a confining harmonic trap with oscillator frequency 
wtr = (27r)100MHz. In Appendix lA2l we provide the min- 
imal Hamiltonian which enables us to describe the scattering 
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FIG. 5: Lowest positive energy adiabatic levels of an interacting 
atom pair near the photoassociation resonance, for the hypothetical 
case of a stable closed channel bound state 7 = 0, as a function 
of the laser detuning A. The light intensity is / = IkWcm"^, 
implying a shift Ao = (27r) 164 MHz. The atoms are asymptoti- 
cally confined in a spherical harmonic trap with oscillator frequency 
wtr = (27r)100MHz. 



properties in the presence of the trapping potential. 

Fig. |5] shows that the lowest positive energy 'scattering' 
state is raised in energy to avoid a degeneracy with the bound 
state approaching from below. Similarly, the states with 
higher positive energy shift upwards such that the adiabatic 
upward ramp of the closed channel state manifests itself as an 
upward propagating 'wavelet' of higher energy level density 

This level structure is modified when the closed channel 
state has a nonzero width 7. As an illustration consider the 
non-hermitian two-level Hamiltonian 
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with all parameters assumed to be real and positive. At 
an avoided crossing of the eigenenergies of i/2icvch where 
E/i = Ciy, the splitting of the levels is real if 7 > 4w and imag- 
inary otherwise. Therefore it depends on the strength w of the 
coupling whether the real parts of the energies avoid cross- 
ing. For sufficiently weak coupling the levels cross, but the 
eigenenergies are both complex, with different decay widths 

(7/2)[l±v/l-16W7)2]. 

As a result of this the real part of the dressed state energy 
may, in the homogeneous case, intrude the positive energy 
continuum as shown in Fig.|3 Fig.|6lshows the adiabatic de- 
pendence of the real and imaginary parts of the dressed state 
spectrum in a 100 MHz trap for a hypothetical decay width 
7 = (27r) 200 MHz. This illustrates in which way the dressed 
state 'intrudes' the spectrum of 'scattering' states with real 
and positive energy. Whereas, at lower energies, the coupling 
is sufficiently weak to allow for a crossing of the real parts, 
scattering states with higher energy couple stronger and there- 
fore avoid to be crossed. The scattering states have, by virtue 
of their resonant coupling to the closed channel state close to 
the curve crossings, a finite lifetime. However, the range of 
detuning over which the scattering states become unstable, is 
of the order of the distance of the harmonic oscillator levels. 
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FIG. 6: Adiabatic levels i5d(A) of the positive energy dressed 
states, as a function of the laser detuning A, for a light intensity 
/ = 1 kWcm~^, for the same trapping frequency as in Fig.|5| The 
width of the closed channel vibrational state was deliberately cho- 
sen on the order of the trapping frequency: 7 = (27r) 200 MHz. 
While the higher levels show avoided crossings as in the stable case 
in Fig.|5| the lowest trap level avoids crossing the rising dressed state 
only in the imaginary part. In the trap the scattering states become 
unstable in the close vicinity of the avoided crossings. 



such that the scattering states are entirely real in energy in the 
homogeneous case considered in Fig.|4] 

We would finally like to note that in Ref. I62ll the Franck- 
Condon overlaps between the different low energy trap states 
and the closed channel bound levels have been studied in de- 
tail. The results corroborate our findings that the coupling 
strengths w vary for the different level crossings as seen in 
Fig.i 



3. Two channel decomposition 

A decisive aspect for the understanding of the two body 
physics close to a scattering resonance is the composition of 
the energy eigenstates in terms of their background and closed 
channel parts as given by Eq. ( I22> . Whereas the closed chan- 
nel component has the functional form of the resonantly cou- 
pled state |0,y), which we have normalized to unity, the back- 
ground channel component depends on W and the energy _Ed 
and therefore varies with A and /. 
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FIG. 7: Fractional population (A/j — 1)/N"l of tlie open cliannel 
component of the dressed state with energy Re£d closest to reso- 
nance. The population is shown as a function of the detuning divided 
by laser intensity (Ao - A)//, for I =0.1 kW/cm^ (dotted line), 
/ = l.OkW/cm^ (dashed hne), and / = lO.OkW/cm^ (solid line). 
Far away from resonance the dressed state is equivalent to the last 
bound state of the background potential such that (A/^ — 1)/N"l ap- 
proaches 100%. Closer to resonance the closed channel state domi- 
nates the dressed state for not too large intensities I. The population, 
however, reaches again 100% at the resonance, where the dressed 
state goes over in a stable scattering state (cf. Sect. lIID"2l . 



In Figure0we show the population 



Nl-l 



^\'\<t>\') = '-^ 



(26) 



of the background channel component of the dressed state as a 
function of the intensity-renormalized detuning (Aq — A)//, 
for / = 0.1 kW/cm2, / = l.OkW/cm^, and / = lOkW/cm^. 
Afd is given by Eq. (I23> . Compare this to the corresponding 
dependence of the dressed state energy and width shown in 
Fig.lU Far from resonance, A ^ Ao(/), where the dressed 
state goes over into the uppermost background channel bound 
state, the population is near 100%. In the intermediate region, 
the population decreases significantly and the dressed state is 
dominated by its closed channel state. 

Closer to resonance, however, the weight of the back- 
ground channel component rises again abruptly. On reso- 
nance, A = Ao(/) (vertical dotted line), the population is 
near 50%, and rises on the positive energy side to 100% 
where the detuning A is such that the imaginary part of the 
dressed state energy vanishes. For / — O.lkW/cm^ (dot- 
ted curve) we have not extended the plot to the detuning 
(Ao - A)// = (27r) 1.12 • lO^MHz kW^^ cm^ where the 
population has risen to 100%. This examples shows that the 
intrusion of the dressed channel state into the positive energy 
continuum as described in the previous section rises with de- 
creasing intensity. In fact, we found that at / = 0.047 kW/cm^ 
the intrusion reaches 6.2 GHz, and below this intensity it is un- 
limited, meaning that for any positive detuning (Ao — A)//, 
a discrete dressed state with an energy with positive real part 
exists, see Fig.|8] 

For large positive (Aq — A)// the imaginary part then also 
approaches j/2. This reflects the complete decoupling of the 
dressed state from the scattering states as one would expect 
from considering the corresponding Franck-Condon overlaps. 




/ [kW cm 



FIG. 8: "Intrusion" of the uppermost dressed state into the posi- 
tive energy scattering continuum. Shown is the maximum energy 
Sd.max ~ R.e(i5d(A, /)) as a function of /, for an appropriately 
chosen detuning A. Below / — 0.047 kW/cm^ the intrusion is infi- 
nite. For I > 0.5 kW/cm^ the maximum intrusion decreases to zero 
with a negative power of I. 
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FIG. 9: Population (A/d — 1)/N2 of the open channel component of 
the dressed state wave function as a function of the intensity / of the 
laser which is tuned to the corresponding (/-dependent) resonance 
frequency, A = Ao(/). For low intensities this population is small 
implying the dressed state to be similar to the closed channel bound 
state. In contrast, strong laser coupling prohibits the dressed state to 
further intrude into the scattering continuum such that it rather re- 
sembles the zero energy scattering state whose closed channel con- 
tribution is negligible. 



In the experiment by McKenzie et al, Ref. |31, flie detun- 
ing of the laser was fixed whereas the intensity I was varied 
in time. Let us therefore have a look at the dependence of 
the background channel population and the dressed state en- 
ergy for the case of different intensities /, while A — Aq (/) 
chosen on resonance. Fig. |5] shows the dependence of the 
on-resonance open channel population on /, as determined 
from Eq. M6\ . At low energies the dressed state resembles 
the closed channel bound state |0^). The higher the intensity 
/, however, the less significant is the closed channel part of 
the dressed state. This may have been anticipated from the 
schematic dressed state spectrum in a trap shown in Fig. |6l 
For sufficiently strong coupling, i.e. large /, the energy RcE'd 
avoids crossing with the lowest scattering state which implies 
that the dressed state goes over into this scattering level as 
soon as it crosses threshold at A = Aq. 
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FIG. 10: Real and imaginary parts of the energy -Ed (A = Ao(/)) 
of the near threshold dressed state on resonance as a function of the 
laser intensity /. As anticipated from Figs.|6|and|9|the dressed state 
is pushed out of the scattering continuum for strong couplings, I ^ 
IkW/cm^ 
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FIG. 1 1 : Open channel component of the radial dressed state density 
for / = l.OkW/cm^ (dashed lines) and / = lOkW/cm^ (solid lines), 
on the PA resonance, A = Ao(/). r and the density are given on 
logarithmic scales. The thin lines show the respective long range 
approximations according to Eq. <27> . The spatial extension of the 
on resonance wave function rises approximately oc I^. 



In Fig.[ro|we explicitly show the dependence of the real and 
imaginary parts of the on resonance energy i?d(A = Ao(/)) 
on the intensity /. The intrusion of the on resonance dressed 
state energy into the positive energy continuum rises with / 
before the dressed state is pushed out again at the highest in- 
tensities, where the situation resembles the case of a stable 
(7 — 0) closed channel state, i.e. Ed -^ 0. 



4. Universal shape 



It was predicted in Ref. 11211 that close to a Feshbach res- 
onance the uppermost bound state is dominated by its back- 
ground channel component which is spatially extended over a 
very large range of the order of the scattering length a. The ex- 
tremely weekly bound molecules formed on the side with pos- 
itive a have binding lengths of unprecedented size, ranging up 
to several thousand OBohr- The production of such molecules 
with ultralong binding length through photoassociation of a 
trapped atom pair has been studied in Ref. L6Z1 . Generalizing 
the arguments in Ref. 14511 to the case where 7 > one finds 
that on interatomic distances r large compared to the van der 
Waals length the background channel component depends on 
r as 
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-r/a{I,A} 



(27) 



From this universal asymptotic behaviour a remarkable fea- 
ture of PA to an unstable molecular state may be in- 
ferred which contrasts the case of a Feshbach resonance 
where spontaneous decay is suppressed. Eq. ([3 yields 
Re[l/a(/,Ao(/))] = abg/|a(/, Ao(/))|2, i.e., accord- 
ing to Eq. (I27> . the spatial extension ((/"d^kl^d^) ~ 
l/(2Re[l/a(/, A)]) of the resonant dressed state is enhanced 
by a factor of |a|/abg — 2S{I)/'y. Consequently, for / = 
10 kW cm~^ the dressed state has an extension which is ap- 
proximately two orders of magnitude larger than for I — I kW 
cm^^. 



Fig. 111! shows the radial probability distribution |(/'(j^(?')P 
associated with background component of the highest vibra- 
tional dressed state. The wave function was determined using 
the minimal Hamiltonian introduced in Appendix |X] for in- 
tensities / = 1.0 kW cm^^ and / = 10 kW cm^'^ of the laser, 
with its frequency tuned on resonance, A = Ao(/). For com- 
parison, also the asymptotic form of the density as approxi- 
mated by the analytic expression in Eq. ( I27> is given. This has 
to be compared to the radial extension of the wave functions 
away, though close to resonance, which is shown in Fig. 4 of 
Ref. l5lll which, at high intensities does not show the extra 
amplification factor. 



III. MANY BODY DYNAMICS OF A BOSE-EINSTEIN 

CONDENSATE UNDER SINGLE COLOR 

PHOTOASSOCIATION 



In the vicinity of a photoassociative resonance the low en- 
ergy scattering properties are modified due to the presence of a 
dressed state near threshold. This has been discussed in detail 
in Sect. HI] For sufficiently high light intensities the modulus 
of the scattering length is significantly increased near reso- 
nance, with in general strong consequences on the dynamics 
of a Bose-Einstein condensed many body ensemble. 

In this section we shall investigate the many body dynamics 
of a trapped Bose-Einstein condensate which is exposed to a 
laser pulse tuned to a fixed frequency near resonance. We will 
focus on possible limitations of the maximum rate at which 
condensate atoms may be associated to diatomic molecules. 
These limits concern, on the one hand, the rates at which 
atoms are lost from the condensate, and, on the other hand, 
the maximum efficiency with which these atoms may be asso- 
ciated to ground state molecules. We will show in detail how 
differences may arise between these limits. 

We will first adapt the dynamical many body theory intro- 
duced in Refs. l ilCl fill |45, 48] to the case of single color 
photoassociation to an unstable closed channel bound state. 
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Central results have been presented in Ref. f51]. 

We shall study the dynamical evolution of the condensate 
as well as of the total number of atoms within the subspace 
defined by their relative motional state being of the form de- 
scribed in Sect. Ill Cl The many body theory allows to system- 
atically treat the condensate mean field and the higher order 
correlation functions. To determine the mean field evolution 
we will take into account the dynamics of the pair correlation 
function. The results of Sect. |ll] allow us to rigorously treat 
the two body physics within the many body formalism. To 
the order of approximation chosen here the many body theory 
includes processes where colliding pairs of condensate atoms 
are transferred to the two body dressed state as well as into 
scattering states with higher relative energy. Moreover, atom 
pairs in the dressed state may be lost from the overall system 
through spontaneous decay caused by the nonzero width 7 of 
the closed channel state. 



A. Microscopic Dynamical Many Body Tlieory 

The many body theory we will adapt for our purposes has 
been introduced in Refs. llOllllll . Its applications presented 
so far include the production of correlated atom pairs in four 
wave mixing of condensates 1 10], the second order correc- 
tion to the condensate ground state energy including the ef- 
fects of three-body collisions 1.5 0.1 . the atom-molecule coher- 
ence in a Ramsey type interferometer lllll . and the condensate 
and molecular dynamics in Feshbach crossing experiments 
145. 481 . A detailed discussion of its application to Feshbach 
resonances within a two-channel picture and of its relation to 
other extended mean field descriptions I I52L 153115411 of the near 
resonant many body physics has been presented in 14511 . 

The theory is based on a resummation of the dynamical 



equations for the single time many body correlation functions 
in terms of cumulants which may be characterized as con- 
nected Green's functions flU BcL 16311 . As the gene ral tech- 
nique have been described in detail before [KJ, ll.^o], we 
shall concentrate on those aspects specific for its application 
to single color photoassociation and required for the subse- 
quent discussion. 



1. Two channel dynamic equations 

The general multichannel many body equations |451 sum- 
marized in Appendix|Clallow to describe a variety of different 
photoassociation schemes in a gas consisting of an arbitrary 
number of different bosonic species. In the following we will 
specify the equations to PA with a single laser the two body 
description of which has been given in Sect.|ll] We may thus 
restrict the many body formalism to the background scattering 
channel coupled via the laser to a closed channel molecular 
state as described in Sect. Ill Al 

In accordance with this picture an atomic condensate mean 
field only exists for ground state atoms. We denote it by 
^(x, t). Since the laser transition is far detuned from the ex- 
cimer potential asymptote, the channel is closed, and a macro- 
scopic population of the excited state may be neglected. To 
take into account the two body correlations which are, in the 
many body system, a measure for the population of the binary 
dressed states lllll . we can restrict the matrix of pair functions 
(cf. Apd.O to its elements corresponding to the background 
channel spectrum, <J'bg(x, y, t), and the closed channel state, 
$ci(x,y,f). 

The mean field for ground state atoms and the pair functions 
then obey the non linear dynamic equations (cf. Eqs. (jCSjl, 
d): 






dt 



dt \ *ci(x,y,i) 



dyvl'*(y,t)(M.g(|x-y|) 



$bg(x,y,t) + 'J'(x,i)*(y,t) +|r!o(t)P.$ci(x,y,t)j, (28) 

(29) 
I 



Here -ffg^(x) = -f?^"^ jlm + Vtrap(x) is the Hamiltonian 
of a single atom in the electronic ground state. Furthermore, 
i7^^(x, y,t) is the general two-channel two-body Hamilto- 
nian of interacting atoms which in the interaction picture and 
rotating wave approximation has the form (cf. Eqs. Q, (|5}): 



i/2B(t) = 
V^2B(t) = 



Tkin 

-?i(A(i)+z7/2)P^ 
P.l^oit) 



+ V2Bit) (30) 
(31) 



H2B may, in general, depend on time through the detuning 

A{t) = u!{t) — LUi, and the intensity of the light I{t) which 
enters the Rabi frequency rio(i). 



The transformation to the interaction picture replaces the 
Hamiltonian in Eqs. ©, (|5} by the Hamiltonian in Eq. ( I30> . 
Thereby, it reduces the coupling VF(R, t) = rJo(0 cos(kR — 
Lot), with R = (x + y)/2, to the envelope J7o(^)/2 and in- 
troduces a transformed closed channel pair function also de- 
noted as $ci- As we chose the zero of the energy to coincide 
with the background potential asymptote, the functions ^ and 
$bg are invariant under the transformation to the interaction 
picture. Physically the transformation of $ci keeps track of 
the fact that the atom pair in the closed channel bound state 
takes up energy and momentum from the photon. Note also 
that we may neglect the phase difference of the light between 
positions x and y since the coupling WP^, acts only over dis- 
tances of the order of the binding length of the closed channel 
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bound state which is much shorter than the wave length of the 
laser. 

For solving the system of dynamic equations ( I28> 
and ( I29> numerically we found it convenient to formally 
solve the inhomogeneous linear equation ( I29> in terms 
of the two body time evolution operator U2B{t,to) = 



Texp{— (i/ft) Jj dt'H2B{t')} and to insert this solution into 
the non-linear equation j28> . The dynamics of the atomic 
mean field is then determined by the single non-Markovian 
non-linear equation: 



ift— *(x,t) 
at 



i/iB(x)*(x,f) 



**(x,i) / dT-^^{K,T) — h{t,T) + {2TThf^^O,hg\V2B{t)U2B{t,to)Mto)) , (32) 

I 



WhCK h{t,T) = i2TThf{0,hg\V2Bit)U2Bit,T)\0,hg)9{t 



t). V2B is given by Eq. ( 13 U . and 9{t — t) is the step function 
that evaluates to zero except for t > r. In deriving Eq. ( I32> 
we have taken into account that the spatial range of the back- 
ground and closed channel potentials is much smaller than the 
scale of the variation of the cumulants. This allows us to elim- 
inate the spatial integrations and leads to the zero momentum 
plane wave state |0, bg) of the relative motion of two atoms 
in the background channel I IL 451 . In the second term in 
the square brackets |$(io)) represents pair correlations at the 
initial time io which will be discussed in following section. 



2. Initial correlations 

The solution of the mean field dynamic equation ( I32> de- 
pends on the initial condensate wave function ^(x, to) as well 



as on the initial pair correlations [^(to))- These are denoted 
as a 'ket' in the sense of (x, x'|$(t)) = |bg)<I)bg(x, x', t) + 
|cl)'l>ci(x, x', t). In general these initial conditions are given 
by the experimental setup. 



For a weakly interacting dilute Bose-Einstein condensed 
gas the non-Markovian mean field equation should recover, 
in the absence of the photoassociation laser, the stationary so- 
lution of the Gross-Pitaevskii equation (GPE). One may fulfill 
this requirement by an appropriate choice of the initial pair 
correlations |$o)- Neglecting the decoupled closed channel 
component $ch we choose the ansatv 



'^b^{Il,r,to) ^ -(2TThfno{Il) 






OT 



(33) 



Here R = (x + y ) /2 and r = x — y are the center of mass and 
relative coordinates, respectively, no is the equilibrium con- 
densate density, and the chemical potential /io is determined 
self consistently from the stationary non-linear Schrodinger 
equation: 



1/2,' \ 



[Tkin + Ftrap(x) + 5(Mo)"o(x)]ny'(x) (34) 



In this equation, the energy dependent coupling constant g in- 
volves the T-matrix associated with the scattering in the back- 
ground channel potential: .g(/io) = (27rfi,)'^Rc[(0|rbg(2/io + 
iO)|0)](cf.,e.g., Ref. |10]).' 

In the dilute gas limit which one has if for the peak den- 
sity the condition ripcakObs ^ 1) holds, Eq. (|33 recovers 
the stationary GPE when g is evaluated at /^o ~^ 0. Note 



that Eq. (I34> . in the limit i ^ oo, exactly follows from 
the dynamic equation (I32> if one inserts the stationary ansatz 



*(x,t) 



1/2 



(x) exp[— i/io(t — to)/h] for the mean field for 



t —^ CX3. 

We finally remark that instead of starting with the ansatz in 
Eq. J33> one could have implemented the stationarity condi- 
tion by choosing the anscitze 



*(R,t) = v/^^;^e-'^« (*-*«)/'*, 
$bg(R, r, t) = no(R) e-*2A'« (*-*«)/'*, 



(35) 
(36) 



for the mean field and background channel pair function and 
derive the initial pair function $bg(R-, r, to) from the formal 
solution of Eq. (|^, cf. Ref. t48J. 
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We underline that the stationarity of the solution for the 
mean field requires, in our dynamical formulation, a care- 
fully balanced equilibrium between the condensate and the 
fraction of the gas described by the pair function. When not 
taking into account the pair correlations in the initial state 
the dynamic equation ( I32> yields, even in the absence of the 
laser coupling, a steadily shrinking condensate fraction. This 
loss rather quickly exceeds the stationary depletion resulting 
from Eq. (I34t . roughly after a time t^ = fi,/|lm/i|. Here 
Im/i = {4:TTh^ahgn/ni){8TTna'l y'"^ — lm[— abg(l — ^fcobg)] 
is the imaginary part of (0, bg|T2s(2/i + *0)|0, bg) evaluated 
to leading order in the momentum k = ^JTrnREfi/h which 
corresponds to the energy of the condensate atoms. 

In the photoassociation experiment 1 8] the peak condensate 
density is of the order of 4 x 10^^ cm^^^, such that the gas, 
without the laser being switched on, is weakly interacting: 
n-peakO-\,„ — 1 ■ 10^^. However, the density is sufficiently 
high that the time scale t^ is on the order of t^ ~ 1.6 ms, 
which is to be compared with the typical duration of 100 /Lts 
of the applied laser pulse. Taking into account the initial cor- 
relations |$(to)) therefore ensures that at low light intensities 
where the induced condensate loss is weak, Eq. J32t yields the 
correct time evolution. 



3. Total density 



In deriving the dynamic equation \32\ for the mean field 
^(x, t) we neglected, in SectionslCland irTlA II the time evo- 



lution of the density F^a of non-condensed atoms, Eq. JC7t . 
In Refs. 1 10, 11] it was shown that this is justified if the ther- 
mal fraction of the gas is negligible at the initial time Iq. 
Moreover, as mentioned in Sect. |Cj in the leading order ap- 
proach a non-condensed fraction arises from the buildup of 
pair correlations, and, in the background and closed channels, 
is given by 



rbg(x,y,i) 
rci(x,y,t) 



dx' 



$bg(x,x',i)$£g(y,x',t) 
$ci(x,x',t)$*i(y,x',t) 



(37) 



The total density of the gas thus is (cf. Eq. ¥2\2V ): 



ntot (x, t) = Tbg (x, x, t) + Tel (x, X, t) + I * (x, t) 1 2 (38) 



Hence, atoms lost from the condensate are transferred to the 
non-condensate fraction in the form of pair correlations $ 7^ 
0. Due to the finite life time of the closed channel state they 
may then be lost from the overall system leading to a decrease 
of A'tot (t) ~ J dx ritot (x, t) with time t. As is shown in 
Appendix |D] the time evolution of the total number of atoms 
is determined by 



l^-W 



dR 



^to)\Ulsit,toM.,c\) ~ {27rhf^^ 



to 



dT^^{R,T)-^e"^'-*h 

OT 



{4>.,c\\ul^{t,T)\0M 



(39) 



where (p{t) — J^ dr (A(r) + 17/2). Note that the term 
including |$(to)) depends on R according to Eq. ( I33> . For 
7 = the total number A^tot would be conserved in time. 



B. Numerical results for photoassociation in ^^Na 



The microscopic theory summarized in Sect. 1111 Al sets us 
in the position to study the many body dynamics of an ultra- 
cold gas under the influence of a photoassociation laser Start- 
ing from a dilute Bose-Einstein condensate of a trapped gas 
of sodium atoms we shall now present the numerical solu- 
tions of the dynamical equations ( I32> and ( I39> . We will first 
consider the condensate evolution and study the effects of the 
non-Markovian character of the interaction term in Eq. (|32} 
thereby concentrating on the achievable photoassociation loss 
rates. 



1. Atomic condensate fraction 



Solving the non-linear non-Markovian mean field equation 
( I32> requires the knowledge of the time evolution operator 
U2B{t, t) which contains all necessary information about the 
two body interactions under the influence of the laser pulse. 
The fact that the equations involve only expectation values of 
this operator and its product with the potential matrix V2B in 
the plane wave state |0, bg) reduces the numerical effort sig- 
nificantly. We shall use the minimal model Hamiltonian in- 
troduced in Appendix IA 1 I to determine the required evolution 
kernels as described in Appendix IA 31 

As initial state we chose a condensate of A^c(^o) = 4.0 • 10^ 
^^Na atoms in a spherical harmonic trap with oscillator fre- 
quency t'ho = 198 Hz. This corresponds to the situation in 
the experiment in Ref. |8|, with the anisotropic harmonic trap 
replaced by a spherically symmetric one with mean frequency 
i^ho = {vxi^yVzY^'^- Solving the GPE gave the initial mean 
field ^(R, to) describing this condensate. The initial back- 
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FIG. 12: (a) Time evolution of the remaining condensate frac- 
tion of the initial atom number A'^(O) =4.0- 10® during a 100 /is 
pulse, with peak intensities of 1.2 kW/cm^ (solid lines), 10 kW/cm^ 
(dashed lines), and 50kW/cm^ (dash-dotted lines). The fraction 
N(t)/N{0) is given on a logarithmic scale. Shown are the solutions 
of the non-Markovian non-linear Schrodinger equation (NMNLSE) 
<32> as thick lines and, for comparison, the evolution according 
to the Gross-Pitaevskii equation (GPE) as thin lines, n(x, t) — 
n(x,0)[l + A'on(x,0)(t - trisc)]~\ triso = 0.5 ^s, with A'o = 
{8TTh/m)\^a{A,I)\. (b) Time evolution of the temporally local de- 
cay 'rate' K{t) — —N{t)/ J dx n (x, i) corresponding to the cases 
considered in (a). The GPE rates are indicated as horizontal lines, 
except for the rate for J = 50 kW/cm^, which is 208 ■ 10"^° cm^/s. 
Time is shown on a logarithmic scale. Since due to the light induced 
shift Ao(/) the resonance condition is first met at the end of the ini- 
tial ramp, t — 0.5 /is, the NMNLSE rates are close to zero before. 
The hatched area indicates the values of K{t) permitted according to 
the limit found in Refs. HHilll, cf. Sect. llIlB2l 



ground channel pair function <I>bg(R, r, to) is then given by 
Eq. ( l33l , with the asymptotic density 7io(R) — |'I'(R, to)P- 

Starting at time to — 0, the laser light, with fixed frequency 
Lu = KeEi^/h+ A, is linearly ramped within 0.5 /is from zero 
to maximum intensity /max, at which it is left for times on 
the order of 100 /is before being ramped back. Corresponding 
to the experimental measurements we determined the fraction 
of atoms [Nc{t) — Nc{to)]/Nc{to) remaining in the conden- 
sate at time t during the laser pulse. The time evolution of 
this fraction is shown in Figure [21 We found that the final 
downward ramp of the laser power has no significant influ- 
ence on the remaining condensate fraction. Fig. E] shows 
the remaining fraction obtained from the non-Markovian non- 
linear Schrodinger equation ( I32t (NMNLSE), for different 



maximum intensities /max, and the laser tuned to resonance, 
A = Ao(/max) = -164MHz/(kW/cm2)x/max- The inten- 
sity /max = 1.2kW/cm^ corresponds to the maximum inten- 
sity reported in Ref. [8]. At this intensity the condensate evo- 
lution can not be distinguished from the time dependence 



nc(x,t) = 



nc(x, 0) 



l+/^o(/max,A)ne(x,0)i' 



(40) 



corresponding to the solution of the GPE with a complex 
scattering length a(/max, Ao(/max)), Eq. (E), which gives a 

loss rate /Co (/max, Ao) = (87rft/TO)|5a(/max, Ao(/max))|- 

Hence, the analysis in Ref. LSJ of their experimental data, in 
terms of the Gross-Pitaevskii theory, is consistent with our 
findings in the framework of the full many body dynamical 
theory. This justifies, a posteriori, to fit our minimal model 
Hamiltonian to experimentally determined quantities as dis- 
cussed in AppendixIXl 



2. Rate limits for condensate loss 



Stimulated by recent theoretical work ISST experimental 
studies of photoassociation have focused on maximum rates 
at which, in principle, condensate atoms may be associated to 
molecules. PA of Bose condensed atoms was, in the past, de- 
tected as loss from the atomic condensate fraction of the gas. 

The upper limit to the possible association rate is set by the 
unitarity constraint to two body scattering |36], which may 
be expressed as the condition that the modulus squared of the 
inelastic scattering matrix element may not exceed one. Ac- 
cording to Ref. 1 8], the density and temperatures in the exper- 
iment were, however, such, that this limit had by far not been 
reached. 

Based on a model of coupled mean fields for the atomic and 
molecular fractions a different limit to the condensate loss rate 
was proposed in Refs. 13a. 14111 . The theory presented in this 
work includes a process through which the molecules, which 
are identified with the closed channel bound states, dissociate 
into the continuum of background channel scattering states. In 
aBEC this limit is of the order of /Cj ^ {h/m)n~^/^, i.e. pio- 
portional to the mean distance between atoms. For the exper- 
iment in Ref. |8] /Cj is far below the unitarity limit, and the 
intensity at which the loss rate is predicted to saturate due to 
this limit |8] coincides with the maximum intensity reported, 
of about /= 1.2kW/cm2. 

We have investigated a range of intensities above / = 
1.2kW/cm2, up to / = SOkW/cm^. Fig.[T2ji shows the time 
evolution of the condensate fraction for on-resonance pho- 
toassociation pulses with /max = lOkW/cm-^ and /max = 
50kW/cm^, together with the coiTesponding time evolution 
according to the GPE. 

We find that, according to our theory which fully takes 
into account the coupling into the continuum of excited states, 
the loss of condensate atoms does not saturate directly above 
/max ^ IkW/cm-^. Nevertheless, the remaining fraction at 
the final time of the pulse is significantly reduced compared 
to the loss derived from the GPE. 
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Already at short times, the loss curves significantly deviate 
from the GPE result. These deviations are due to the non- 
Markovian character of the evolution equation ( I32> . The on- 
resonance scattering length Q is, for high intensities, domi- 
nated by its imaginary part Ima(/, Ao(/)) = — 2abg(S(/)/7. 
The Markov time of the many body evolution is of the order 
of the duration of a two body collision which is comparable to 
m\a\'^ /h. For high light intensities, where \a\ becomes large, 
this timescale is significantly larger than the time over which 
the condensate mean field changes. 

In Fig. [21 the non-Markovian character of the mean field 
evolution emerges as a temporally varying decay 'rate' 



K{t) = -N,{t) 



dx.n^{x, t) 



(41) 



In the Gross-Pitaevskii theory the nonlinear interaction term 
is local in time such that K is a. constant, as indicated by the 
horizontal lines in Fig.^] As soon as the interactions are res- 
onantly enhanced, at around t = 0.5 /is, where the intensity 
dependent shift Ag (/) starts to compensate the laser detuning 
A, the rate starts to rise. Initially, however, K{t) rises itself 
at a maximum rate (note the logarithmic time scale). Whereas 
for /max = 1.2kW/cm^ the rate K{t) quickly saturates at 
the Gross-Pitaevskii value Kq, for higher intensity the non- 
Markovian time evolution induces larger fluctuations before 
K{t) settles in at a final value. In view of these oscillations 
it is not justified to characterize the decrease by a single loss 
rate. In Sect. lIIIB4l we shall discuss in more detail the physi- 
cal grounds for these fluctuations. 

Fig. 1 12| allows to relate our results to the limit proposed 
in Refs. f3»,M7|. Inserting the local density nc(R, i), calcu- 
lated with Eq. (|32|l, into X,i(R,t) = {h/m)[nc{Tl,t)]-~'^/^ 
we calculated the limit K{t) given by Eq. ( 14 H . with 
Nc{t) — — J (iRA'j(R, t)nc(R,t), for a maximum intensity 
of /max = 50kW/cm^. The values of /ir(t) allowed according 
to this limit are indicated by the hatched area in Fig. [21 Even 
for such high intensities for which the condensate density de- 
creases rapidly, thereby raising the upper limit Kj, we find 
that the loss determined by Eq. (j^D considerably exceeds this 
limit. 

In summary: At the light intensities applied in the experi- 
ment in Ref. [8] the loss of atoms from the condensate may be 
described, to a good approximation, by the GPE with a com- 
plex coupling constant proportional to the scattering length 
given in Eq. Q. At higher intensities, the dynamics may no 
longer be described by a single rate constant. Moreover, the 
fraction of atoms lost does not increase as strongly as pre- 
dicted by the Gross-Pitaevskii theory. However, the limitation 
of the loss rate is, according to our results, not as stringent as 
predicted in Refs. Mlilll . 



3. Dissipative molecule formation 

The results of the previous section have shown that in the 
vicinity of a photoassociative scattering resonance the con- 
densate dynamics may deviate significantly from basic mean 
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FIG. 13: Time evolution of the fractions Ni{t)/N{Q) of atoms re- 
maining in the condensate (i = 'c') and in total {i = 'tot'), of the ini- 
tial atom number A'^(0)= 4.0 ■ 10® during a 100 /is pulse, with peak 
intensities of 1.2kW/cm^ (solid lines), lOkW/cm^ (dashed lines), 
and 50kW/cm^ (dotted lines). The fractions Ni{t)/N{0) are given 
on a logarithmic scale. The respective higher lying lines correspond 
to the total atom number, the lower ones to the condensate fraction. 
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FIG. 14: Loss fractions 1 - A'^i(tfln)/A'^c(0), i = 'c', 'tot', of the 
condensate (filled squares) and total atom (filled diamonds) loss at 
the end of the laser pulse, ifln = 100 /is, for different laser intensi- 
ties /max. The fraction indicated by the diamonds corresponds to the 
number of ground state molecules formed. The solid line shows the 
result of a solution of the GPE in local density approximation. The 
inset shows the same quantities on a logarithmic scale. The modifica- 
tions of the curves compared to Fig. 3 of Ref. 1 5 1 ] arise from taking 
into account, here, initial correlations as described in Sect. lIIIA2l 



field behaviour, in particular at light intensities where the 
modulus of the scattering length Q is significantly enhanced. 
In Section IIIDI we considered in detail the properties of the 
discrete near resonant dressed state with the highest energy 
Re/^d- We shall now show how properties of this state mani- 
fest themselves in the many body dynamics of near resonant 
photoassociation. An observable of interest in this respect is 
the time evolution of the total number Ntot{t)- We recall that 
this number may decrease due to the finite life time 7~^ of the 
closed channel state. 

Using the solutions of Eq. ({32} for different maximum laser 
intensities /max and detunings A = Ao(/max), as well as 
Eq. ( I33t and the two body time evolution U2B determined 



with the minimal model Hamiltonian given in Appendix 1X1 
we calculated Ntot{t) fromEq. ( I39> . Fig.ll3lshows the frac- 
tion 1 — iVtot(i)/-^tot(0) of atoms remaining at time t of the 
initial total number of iVtot (0) = Nc{0) = 4.0 • 10^ sodium 
atoms, for the same light intensities as in Fig. [2] For com- 
parison we add the condensate time evolution from Fig. El 
which corresponds to the respective lower lines. 

Fig. nH shows the final fraction 1 - iVtot(ifin)/^tot(0) of 
atoms remaining at t — tun — 100 fis (filled diamonds, lin- 
early interpolated by a dashed line). Again, we compare this 
fraction to the final condensate fraction 1 — A'c(ifin)/^c(0). 
At low intensities the two fractions are identical, i.e. all atoms 
missing from the condensate have been lost through the reso- 
nantly coupled unstable closed channel molecular state. How- 
ever, when increasing the intensity above ca. 1 kW/cm^ the to- 
tal loss of atoms saturates and decreases slightly at the highest 
intensities. According to Eqs. (I37> and ( I38> the loss of con- 
densate atoms is due to the formation of binary correlations in 
the gas. 

It was shown in Refs. llllll45ll . for the formally similar case 
of a Feshbach resonance, that the pair functions $bg, "I'd con- 
tain the information about the population of the background 
and closed channel components of the two body bound and 
scattering states. With the corresponding quantum mechani- 
cal wave functions serving as form factors the pair functions 
yield observables for the number of molecules and excited 
atom pairs within the many body emsemble. These observ- 
ables were measured in a number of recent experiments where 
molecules and excited or 'burst' atoms are produced using 
Feshbach resonances ITll 111 145-48.1 . 

In contrast to Feshbach resonances, a meaningful observ- 
able measuring the number of molecules produced in single 
color PA is neither given by the number of atom pairs in the 
dressed state with the highest energy Re-Ed which has been 
discussed in Sect. IIIDI nor by its closed channel component 
(cf. Sect. IIID3> . Due to the fact that the decay time of the 
closed channel state is much shorter than any other timescale 
near resonance, any molecules formed are bound states of the 
background channel ground state potential. Since we do not 
explicitly keep track of the target states of this decay pro- 
cess a meaningful observable for the maximum number of 
such molecules formed is the decrease of the total number 
of atoms A^tot- The diamonds in Fig. ^Jmay therefore be 
seen as marking an upper limit to the number of molecules 
formed from the condensate atoms, which saturates at inten- 
sities above /,„ax ^ IkW/cm^. 

We compare Fig. O with the intensity dependence of the 
discrete dressed state close to threshold, shown in Figs.|5]and 
ITOl As expected, the decreasing significance of the closed 
channel component of the dressed state above ca. /max — 
1 kW/cm^ results in a reduction of the efficiency with which 
the condensate atoms may be associated to ground state 
molecules. The two body correlations induced by the high 
intensity laser partly relate to dressed two body states with 
a strongly enhanced lifetime. According to Fig. [TO] the life- 
time of the dressed state at /max = lOkW/cm^ is on the order 
of 16 /is, a factor of 2000 larger than the lifetime of the bare 
closed channel vibrational state \6,,). 
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FIG. 15: Time evolution of the radial density profile 

7ii(r, t)/nc(0, 0) of the fraction of condensate atoms (i —c) and the 
total number (i =tot) remaining at times f = 1, 10, 100 /is, for three 
different intensities / — 1.2 (solid), 10 (dashed), and 50kW/cm^ 
(dotted lines). At time i = all evolutions start with a conden- 
sate of 4 ■ 10® sodium atoms in a spherically symmetric, harmonic 
198 Hz trap. The respective upper curves correspond to the total 
atomic density distribution ntot , the lower curves to the condensate 
profile Tic- All profiles are normalized to the initial peak density 
ripcak ~ nc{0, 0). Note the logarithmic scale in the lower diagrams. 



While the loss of atoms by spontaneous emission is reduced 
at high intensities, the laser stimulated loss to excited non- 
condensed atom pairs is gaining weight. This explains why 
the saturation in the total atom loss is not accompanied by 
an equivalent saturation in the condensate loss. Despite the 
fact that the closed channel component of the near threshold 
dressed state becomes negligible the resonant coupling to this 
state still allows for the population of the scattering states of 
atoms asymptotically in the internal ground state. This trans- 
fer to the non-condensed fraction of the gas is a characteristic 
result of the full dynamical many body theory described in 
this article. 



4. Collective excitations 

We would finally like to present results for the time evo- 
lution of the density profile of the atom cloud trapped in a 
spherically symmetric harmonic potential. For the calcula- 
tions presented in the previous sections we have taken into ac- 
count, in Eq. ( I32> . a trapping potential Vtrap = i^ho^^' ^^'^^ 
Who — (27r)198Hz, corresponding to the mean trapping fre- 
quency in the experiment of McKenzie et al. |8]. As we have 
seen in Figs.ll2h.ll2b the condensate time evolution exhibits 
considerable oscillations at higher light intensities. A detailed 
study of these oscillations is beyond the scope of this paper 
and will be presented elsewhere. The oscillations are, how- 
ever, expected to be understood analogously to a similar effect 
in a Feshbach resonance experiment of the kind presented in 
Ref. t3- It was shown in Refs. OQ that due to the fact 
that the condensate atoms and molecules are not completely 
orthogonal to each other close to a Feshbach resonance, the 
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corresponding fractions of the gas exhibit residual oscilla- 
tions. In the dynamical evolution presented in Refs. 1 ,11. .121 
these oscillations are seen during the evolution period between 
the initial and final Ramsey pulses. Their frequency largely 
corresponds to the binding energy of the molecules formed on 
the positive scattering length side of the resonance. Close to 
resonance, however, the experimental measurement of the os- 
cillation frequency exhibited a shift 1 14| which may be under- 
stood as caused by the many body environment in which the 
molecule association occurs |66]. Our discussion of the near 
resonant dressed state spectrum in Sects. InDTI and HI D2l has 
shown that for the laser detuning A = Ao(/max) during the 
pulse the real part of the uppermost dressed state is positive, 
of the order of a hundred kilohertz for /max = lOkW/cm^. 
The time evolution of the rate shown in Fig. 112b thus indicate 
a many body induced frequency shift which enables oscilla- 
tions despite a positive dressed state energy. 

Fig. E] shows the radial density profile of the conden- 
sate fraction, \'^{R,t)\'^, and of the total number of atoms, 
ntot{Rit), for three different intensities /,nax of the light, at 
various times during the light pulse. The density of both, the 
condensate and the total atom number develops spatial oscilla- 
tions, especially at the highest intensities. These collective ex- 
citations arise because the frequency of the temporal oscilla- 
tions depend on the local density of the cloud, which provides 
a further indication of their many body origin. Our results 
suggest that an experiment with a Ramsey interferometer like 
setup as reported in Ref. lll3ll may be used for precise molec- 
ular spectroscopy and to study the influence of the many body 
environment of the ultracold quantum degenerate gas. 



IV. CONCLUSIONS 

A dynamical theory of photoassociation (R\) in a Bose- 
Einstein condensate has been presented which sets us in a po- 
sition to study situations which due to the external conditions 
are no longer correctly described in the Gross-Pitaevskii ap- 
proach. The single color PA process type which is the subject 
of this article is described in a two-channel picture, account- 
ing for an open background scattering channel and a molecu- 
lar vibrational state in a closed channel potential of interact- 
ing electronically excited atoms. The two-channel description 
was chosen since it allows to include the decay of the closed 
channel state due to spontaneous photon emission and its con- 
sequences for the many body evolution in a straightforward 
and convenient way. We point out, though, that a single chan- 
nel picture in principle suffices to obtain accurate results for 
the many body dynamics of the condensate fraction, as has 
been shown in Refs. 1 11, 481 for the case of a Feshbach reso- 
nance. 

The modifications of the discrete spectrum of binary bound 
states near the photoassociation resonance have been de- 
scribed and the role of the finite liferime of the closed channel 
state has been shed light on. We have extended our minimal 
Hamiltonian description of the low energy two-body scatter- 
ing properties to the case where the atom pair is additionally 
confined by a tight spherically symmetric harmonic trap. As 



an example, the adiabatic level structure near a Feshbach res- 
onance as obtained from this approach was compared to cor- 
responding results in Ref |68] from numerical multichannel 
calculations and an effective scattering length model. Excel- 
lent agreement confirmed the validity of the approach devel- 
oped here. By applying the tight trap Hamiltonian to single 
color PA we could illustrate the level structure behind the near 
resonance dressed state behavior found in the homogeneous 
case. The real part of the dressed state energy was shown to 
intrude the positive real energy continuum of scattering states 
by crossing a number of low energy scattering levels before 
settling in and adiabatically going over into a scattering state 
itself. Thereby the scattering level energies become complex 
close to the crossings, within a window of the order of the trap 
level spacing. In other words, when tuning through the reso- 
nance the uppermost dressed level may survive on the side 
where the real part of its energy becomes positive, thereby 
avoiding a crossing of the scattering levels only in their ener- 
gies' imaginary part. We have discussed, furthermore, univer- 
sal properties of the dressed state closest to the open channel 
threshold and shown that the spatial extent of its wave function 
on resonance scales, for high laser intensities /, quadratically 
with /. This is a property specific to the case of an unstable 
closed channel state, i.e., it is not present in the bound state 
near a Feshbach resonance whose binding length scales with 
the scattering length and therefore linearly with the width of 
the resonance. 

On the basis of the two body description developed here we 
have studied the dynamical evolution of a trapped condensate 
of sodium atoms whose interactions are strongly modified by 
a laser which is rapidly switched to a high light intensity, com- 
parable to the situation in the experiment in Ref. 1 8]. We have 
computed the temporal evolution of the condensate fraction 
as well as of the total atom number for different intensities at 
the verge and above those reached in experiment. The results 
exhibit the significance of the beyond mean field approach. 
A simple loss rate description as in the Gross-Pitaevskii pic- 
ture is not possible at high light intensities, and the amount of 
atoms lost within a certain time interval from the condensate is 
considerably limited compared to what may be expected from 
two-body theory. Hence, our results corroborate earlier pre- 
dictions of such a rate limit, but do not confirm the predicted 
value of the limit by allowing a faster loss. We conclude that 
at the intensities reached in the experiment of McKenzie et 
al. (^ a saturation of the condensate loss could not have been 
observed, but should be found if the intensity could be in- 
creased by roughly an order of magnitude. Similar studies for 
other PA resonances are beyond the scale of this paper but 
may reveal saturation due to many body effects at more real- 
istic light intensities. The predicted condensate loss rate lim- 
itation is due to incoherent production of excited atom pairs, 
of the kind observed as a burst in experiments with Feshbach 
resonances, cf., e.g., Ref. lOl . 

We have also found a saturation of the number of ground 
state molecules which may be formed in the spontaneous elec- 
tronic deexcitation process of the decaying closed channel 
bound states. This saturation may be understood as result- 
ing from the towards high intensities decreasing significance 
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of the closed channel part in the two component dressed state 
vector vs. its open channel part. Similarly as the uppermost 
bound state on the positive scattering length side of a Fesh- 
bach resonance 1 12] the extremely long range dressed state 
lives almost exclusively in the open channel, rendering the 
closed channel population irrelevant when describing the dy- 
namics in a two channel picture. For high light intensities, 
the achievable number of molecules is predicted to be signifi- 
cantly below the number of lost condensate atoms. 

We have finally illustrated the collective excitations formed 
in the atomic cloud by the high intensity laser pulse. These 
arise from a temporal oscillation of the condensate density 
during the laser illumination whose frequency depends on the 
local value of this density. These oscillations, the detailed 
analysis of which is a subject of future work, are expected to 
be caused by a mean field shift of the dressed state energy 
and to be understood similarly as the residual coherent oscil- 
lations between condensate atoms and condensed molecules 
as observed in the Feshbach Ramsey interferometry 1131 and 
described in Ref. lOl . 
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APPENDIX A: A MINIMAL MODEL FOR TWO BODY 
SCATTERING 

The scattering of a pair of atoms at very low kinetic energies, 
as in a Bose-Einstein condensed ensemble, is, to a very good 
accuracy, determined by a few parameters only. As discussed 
in detail in Sects. HTCl a complete set of such parameters fix- 
ing a minimal two channel model Hamiltonian is given by the 
background scattering length Obg, the energy E^i of the last 
vibrational level in the background potential, the width 6 {I) 
and shift Aq (/) of the laser induced resonance, and the decay 
width 7 of the resonantly coupled closed channel state. Any 
Hamiltonian which reproduces these quantities may serve as 
a precise model to describe the full low energy scattering. 

In this appendix we recall the basic ideas of the separable 
model introduced and discussed in detail in Refs. JllL 14511 and 
present the separable parameters used throughout this work to 
calculate both the static and dynamic scattering properties in 
Sects, nilandlllll An extension for the case of scattering in a 
spherical trapping potential follows which will be used in Ap- 
pendix|n]to study the adiabatic energy level structure. Finally 
the method for determining the two body time evolution will 
be described. 



1. Two channel model for asymptotically free atoms 

Separable models are related to the observation that in cer- 
tain scattering problems the transition matrix is to a very good 
approximation determined by a single pole | HJ 16711 . For in- 
stance, in the case of low energy scattering considered here, 
this pole corresponds to the last vibrational state of the poten- 
tial. Knowledge of the potential and therefore of the spectrum 
of bound states would allow to calculate the T-matrix in this 
pole approximation. While such a procedure would require an 
intermediate step with the much larger amount of information 
given by the potential the separable approach may be used in 
the opposite direction to fix an ansatz for the last vibrational 
state of the potential V in terms of a minimal set of defining 
parameters, cf. Ref. 1 11]. Since the detailed form of the po- 
tential is irrelevant for collisions in ultracold gases, the low 
energy scattering properties are fully determined by this set. 

In the single channel description of a scattering resonance 
presented in Ref. LI 1.1 it is sufficient to provide the bound state 
energy and the scattering length for any configuration close to 
and away from the resonance. The results of Ref. \IY\ have 
demonstrated that such a single channel description is suffi- 
cient to reproduce the relevant low energy scattering prop- 
erties and physical phenomena in the vicinity of a Feshbach 
resonance. In Refs. |45, 48] the single channel approach has 
been compared to a two channel description of a Feshbach 
resonance, and their results confirm that a single channel de- 
scription suffices for the applications considered. 

For practical convenience we will use a two channel de- 
scription for single colour photoassociation. The reasons for 
this are the finite life time of the closed channel state and the 
temporal variation of the light intensity and therefore of the 
width and shift of the resonance. Thereby the two channel 
model automatically provides the correct dependence of the 
complex energy of the last vibrational state on detuning and 
intensity of the laser. Our model consists of a separable ansatz 
for the background potential. 



Vbg = |Xbg)^bg(Xbg| 
and the laser coupling 



(Al) 



(A2) 



with real parameters ^bg and C^. Since the resonance shift 
Ao(/) and width 5{I) (cf. Eq.Q) depend linearly on / at the 
intensities considered we have pulled a factor of \/l out of C,. 
The practical convenience of the above ansatz results from 
choosing a Gaussian form for the states: 



(rlxbg) = exp(-rV20/(%/2^abg)3 (A3) 

(r|x> = exp(-rV2cr2)/(\/2^cr)3. (A4) 

Here, the parameters crbg and a define the width of the Gaus- 
sians and thus reflect the spatial extension of the van der Waals 
tail of the background potential and of the closed channel 
molecular state, respectively. 

While the closed channel decay width 7 enters the 
Hamiltonian given in Eqs. @ and (|5} as a real parame- 
ter, the four other parameters fix the separable ansatz of 
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Eqs. JA1> and ( IA2> . For the photoassociation resonance 
shown in Fig. [2 the parameters assume the following val- 
ues: m^hg/{4:TTfi^) = -191.4 aBohr, crbg = 23.97 asohr, 
e = -(27r)57GHz (aLhr/(kW/cm2))i/2, ^ ^ 3.55aBohr, 
and 7 = (27r)18.5MHz. These values have been determined 
as follows. 

Inserting the separable form ( IA1> into the Schrodinger 
equation for the last bound state of the background Hamil- 
tonian yields an analytic expression of E-i in terms of 
fbg and (Tbg: Multiplying the integral form {(p-i) = 
Go{E-i)Vhg\(p-i) of the Schrodinger equation by (xbg| from 
the left gives 



l-ebg(Xbg|Go(i?_i)|Xbg)=0, 



(A5) 



The separable ansatz allows to calculate the expectation val- 
ues of Go with respect to the states |xbg) , |x) ^s 



{Xa\Go{E)\xb) 



m 



4:T:h?(7ab 



^ah^ 



''eiic{xab) 1= 



(A6) 



Here Xa,Xb e {xbg,x}, o-Qb == y/{al + af)/2, (Ja,(Tb e 
{cbgiCr}, and erfc(x) = 1 — 27r^^/^ /^^^ dw exp[— u^] is 
the complimentary error function of the argument Xab — 
V-mEaab/h. 

A second relation between the background scattering 
length, the energy _E_ i , and the separable parameters ^bg and 
tTbg is obtained from the background T-matrix which, in the 
separable potential approach takes the form lllll 



rp ( \ ^ IXbg)Cbg(Xbg| 

"'^ ' l-ag(Xbg|Go(z)|xbg) 



(A7) 



Using the definition iTrfi^abg/m = (27rfi.)3(0|rbg(0)|0) of 
the background scattering length in terms of the background 
T-matrix T'bg(z) = Vbg + VbgGo(2)Vbg one obtains an ex- 
pression of flbg in terms of ^bg and (xbg|Go(0)|xbg)- With 
Eq. (IA6> this gives; 



Cb 



Airh^ahe 



1 



Qbg 
V^CTbg 



(A8) 



This equation yields, inserted in Eq. ( IA5> and using Eq. ( IA6> . 
a relation which determines (Jbg in terms of abg and E^ i : 



for the scattering states in the background channel and the ex- 
pression 



Gbg(z) = Go(z)[l + Tbg(z)Go(z)] 



(All) 



for the background channel Green's function, Eq. dlOt . in 
terms of the background channel T-matrix, Eq. JA7> . one ob- 
tains: 



dAp/dl 
dS/dl 



1 



1 



^/7^a l + (CTbg/o")^ 



V^<T V l+(<Tbg/<T)^ 



2TTh3 as 

ni dl 



1 



"bg / 2 

^ay l + (<Tbg/^)" 



(A12) 



(A13) 



To arrive at these relations we have also used Eq. (IA6> and 
the approximative linearity in / of 5 and Aq. Eqs. (IA12t and 
( IA13l l. together with Eqs. ( IA8> and ( IA9> determine the pa- 
rameters of the minimal model Hamiltonian given above in 
terms of experimentally known or otherwise theoretically de- 
termined quantities. 

In Sect. Ill D1 we used the separable model to calculate the 
near resonant bound states and their energies. As an example 
we give the equation which was used to determine the energy 
dependence of the dressed states on the intensity / and detun- 
ing A. Using Eqs. (|24|i, (|A2li, (|X3, and dATTl we find: 



-Ed = A - i-f/2 + C^I 



(x|Go(i?d)|x> 



ag((x|Go(gd)|Xbg))^ 
l-Cbg(xbg|Go(-Bd)|xbg) 



(A14) 



Here we have also used, that the time independent left eigen- 
state (0^ I may not be distinguished from the corresponding 
right eigenstate {4>u\, such that {(f)u\W = C,"^{x\ (cf- Ap- 
pendix|E}. The expectation values of the Green's function Gq 
may be determined using Eq. ( IA6t . 



2. Extension to trapped atoms 



1 _ ^:^e^^erfc(x) = 0. 

(Tbg 



(A9) 



Having determined (7bg in this way, Eq. ( IA8> finally yields 

^bg. 

Introducing the pole approximation with respect to the 
closed channel state \(j)^) we have already obtained the off 
diagonal coupling in a separable form. The parameters ^ and 
(7 are therefore directly determined from Eqs. ( I19> and ( I20> . 
Using the Lippmann-Schwinger equation 

|0(,+'> = |p> + Go(pVm + iO)Tbg(pVm + iO)|p) (AlO) 



Within an external trapping potential the relative motion of 
two atoms is additionally restricted with the effect that the 
spectrum of positive energy scattering state is no longer con- 
tinous. In this section we introduce a separable potential ap- 
proach to treat low energy two body scattering within a con- 
fining spherically symmetric harmonic trap. For the purpose 
of this article our goal is to calculate the spectrum of nega- 
tive and positive energy dressed states of the relative motion 
within the trap and their dependence on the detuning and in- 
tensity of the laser To this end we will fix the separable pa- 
rameters such that in the homogeneous limit they yield the 
known low energy scattering properties. 
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a. Two channel separable potential 

An effective scattering length model of ultracold atomic col- 
lisions in tight harmonic traps has been discussed in detail in 
Ref. 168]. Like this model our minimal Hamiltonian model 
with separable interactions accurately reproduces the positive 
energy spectrum of the atomic relative motion determined 
from multichannel closed coupling calculations presented in 
ll68ll . We expect the accuracy of our model to hold for tight 
traps as long as the oscillator length considerably exceeds the 
length scale /vdw of the van der Waals interactions. 

We consider two atoms whose relative motion in the back- 
ground channel is, without the coupling laser, described by 
the Hamiltonian 

i^2B=i/o°|bg)(bg| + |(^,,cl)^,(0,,cl|+l^2B. (A15) 

Here, the second term is equivalent to the corresponding term 
in Eq. @, and V2B is the two channel interaction picture po- 
tential matrix given in Eq. (|5|i. The Hamiltonian of trapped, 
non-interacting atom pairs reads, in the spatial representation: 



m°{r 



h^V^ 



1 2 



(A16) 



where u>\^oI2-k is the harmonic oscillator frequency. The mo- 
tion of non-interacting atoms is restricted to the spherical trap 
modes |n, /), with 



a, 



hoi 



r(0)| 



|n,0=</KO, 

En,l — fi^ho 



-+2n + l 



(A17) 
(A18) 



Since the trap is isotropic we may, as in the uniform case, 
consider only the s-wave term (l ~ 0) in the expansion of the 
scattering amplitude. We define an s-wave scattering length 
by 



«bg(^) = ^^AA-i(o,o|rbg(i?)|o,o), 



(A19) 



where |0, 0) is the (n = 0, / = 0) oscillator ground state de- 
fined by Eq. ( IA17t . and the T-matrix is evaluated at zero or 
positive energy E > 0. Furthermore, J\f = {2'k)^-^/^1^^, 
determined by the oscillator length /ho = y/h/nuJhxi, is a 
normalization factor which is chosen such that in the homoge- 
neous limit the scattering length ab°(0) defined by Eq. ( IA19> 
is equivalent to the usual Cbg. 

We introduce a separable potential as in Eqs. (IA1> . (IA7> 
and make the following convenient ansatz for the trap repre- 
sentation of the state |xbg) in terms of a real parameter ahg'- 



("•,0|xbg 



{l + 2n)U 



1/2 



,/iL 



(A20) 



il) =T{a + l)/[r(6 + l)r(a - b)] denotes the generaUzed 
binomial coefficient. 

Using the spectral expansion 



Go{z) 



Z^ 



\n,l)iz-Ell>})-\nJ\ 



(A21) 



of the free energy dependent Green's function its expectation 
values with respect to the states Ixbg), \x) ^rc (cf. Eq. jA6» : 



{Xa\Go{E)\xb) 
m 



47r?i20r/ho 



I -at 



2/- 2 ) J(y^^o.b) 



(A22) 



Here, y 



-E/huJho, a-ab = exp[-cr^j//^J, aab 



^/{al+'af)J2, (Ta.cTb e {(Tbg, o-}, and the integral J(y, a) 
may be expressed in terms of a hypergeometric function: 



J{y, a) — v2 / du 
Jo 

2\/2 



,1/2-y 



\/l — u^a^ 



■i-2y 



F(i,i(|-y);i(|-y);a2). (A23) 



Inserting the separable T-matrix, Eq. (IA7> . into Eq. ( IA19l l 
and using the ansatz ( IA20t we obtain (cf. Eq. ( IA8» : 



ag(i?) = 



A^h?all{E) 



m 



1- 



<1{E) / 2 



'-^A.e-<"h 



(A24) 



One may show that for any fixed — y = E/Hloy^o > 0, 
Eq. JA24> is, in the weak trap limit Who -^ 0, equivalent 
to Eq. ( IA8> . For the following considerations we define ^bg 
through Eq. ( IA19> with E ~ huj\^o/2. 

Using Eq. JA22> we then obtain a relation determining the 
spectrum of energies E in the background channel analogous 
to Eq. (|A9}: 



;S:G-»)''<--""-"^^'-' 



(A25) 



where y = -E/fiwi,o, Obg = al°^{hjJi,o / 2) . Again, in 
the weak trap limit, Who -^ 0, this relation is equivalent to 
Eq. jA9> . We therefore used Eqs. jA8> and jA9> to fix the 
separable parameters ^bg and (Tbg. 

The separable parameters a and (^ (cf. Eq. (IA2> ') are de- 
termined in the same way as in the uniform case. Using the 
arguments which have lead to Eqs. (IA12t . ( IA13> we find: 



(A26) 
(A27) 



d^o/dl ^ |l+ebg(x|Gbg(0)|x)P 
dS/dl 



(xlGbg(O)lx) 
d5/dl 



L(x|Gbg(0)|x) 



1/2 



The expectation values of the Greens function Gbg(O) are 
evaluated using Eqs. JA7l i. ( lAlU . and ( IA22t . 
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b. Accuracy of separable approximation in tight traps 

To check the accuracy of the minimal model introduced here 
we have considered the scattering of ^"^Na atoms in the {F = 
l,mp — +1) stretched state, in the vicinity of the Fesh- 
bach resonance at Bq = 90.1 mT for the case that the atoms 
are confined within a tight spherically symmetric trap with 
Who = (27r)500kHz. In Ref. 1.68.1 the five lowest positive en- 
ergy levels were presented as a function of the magnetic field 
detuning from the Feshbach resonance, determined with mul- 
tichannel close coupling methods from high precision data for 
the relevant interaction potentials. Here, we use Eq. ( I24> . with 
A(B) == {dE^/dB){B - Bq) = (52.4 MHz /i/mT) (5- 
Bq) and 7 = 0, together with the separable ansatz (|A20} 
for the interactions modified by the trap, to determine the 
spectrum of positive energy states. Inserting the parame- 
ters m^bg/(47r/i^) = -184aBohr, crbg = 26.5aBohr, C = 

~44.0GHzft,a3Qjjj. and a = 17.4aBohr which were deter- 
mined from the low energy scattering properties of untrapped 
atoms close to the resonance (cf. Ref. |48|), we obtain the de- 
pendence of the energy levels on B. The results are shown 
in Fig. [^('h-' -signs) where they are compared to the results 
from Ref. |68], of numerical calculations (filled circles) as 
well as of an effective scattering length model (solid line). 

In Ref. 1 69] the energies E of the relative motion of trapped 
atom pairs were found, using the pseudopotential approxima- 
tion to the binary interactions, to be given by the solutions of 
the equation 1701 



4500 



1 



ho 



V2 



tan 



TT /I 
2 V2 



r(i 



(A28) 



where a is the scattering length entering the pseudopotential 
and y = —E/huiho- Extending this approach, the authors of 
Ref. 1 68 ] argue that a needs to be replaced by an effective scat- 
tering length UesiE) = — ta,n[So{k)]/k, where E — fPk'^/m. 
and (So(fc) is the momentum dependent s-wave phase shift. 
This causes the magnetic field value, where a particular energy 
level corresponds to a diverging scattering length, to be shifted 
for subsequent levels, as indicated by the slanted dashed line 
in Fig. ^] For the level structure determined from Eq. ( IA28> 
the dashed line would be essentially vertical. 

As discussed in detail in Refs. Illl 14511 the parameter a^g 
introduces the scale of the van der Waals length Ivdw into 
the description of the atomic interactions. Hence the pseu- 
dopotential approximation should be recovered in the limit 
CTbg -^ 0. Using the relation F{a, b; c; 1) = r(c)r(c — 
a — b)/[T{c — a)r(c — b)] between the hypergeometric func- 
tion and Euler's F function one may indeed show that in this 
limit Eq. ( IA28> is equivalent to Eq. (IA25> . identifying thereby 
a = all. 

Hence, the excellent agreement of our results with the re- 
sults of Ref. |68|, as seen in Fig.[^ not only shows that the 
minimal model Hamiltonian remains applicable outside the 
Wigner law regime, where the scattering length is much larger 
than the harmonic oscillator length. It also nicely illustrates 
that the particular trap induced level shifts close to resonance 
are strongly influenced by the far ranged van der Waals inter- 



4000, 

3500 

3000, 

B 2500 







90 90.02 90.04 90.06 90.08 90.1 90.12 90.14 90.16 90.18 90.2 
B[mT] 



FIG. 16: Comparison of adiabatic energy levels calculated in this 
paper with the result from Ref. 1 68] for the five lowest positive energy 
levels near the Bo = 90.1 mT Feshbach resonance of ^^Na atoms 
colliding in the (F = \, mp = +1) hyperfine ground state, which 
are confined in a i^ho = 500 kHz spherically harmonic trap. I thank 
P. Julienne for providing me with Fig. 4 of Ref. I6m1 . 



actions. 



3. Determination of the time dependent coupling function 

The essential ingredient of the non-Markovian non-linear 
equation ( I32> for the mean field \E'(x, t) is the antiderivative 
h{t,T) = (27r?i)3(O,bg|F2B(i)C/2B(i,r)|O,bg)0(t-r)ofthe 
expectation value of the time dependent T-matrix between the 
zero momentum plane wave states in the background channel. 



dh(t,T)/dt = g{t,T) 



-G2B(t,r)F(r)]|0,bg). 
(A29) 



In the following we will briefly summarize the procedure in- 
volving the minimal model Hamiltonian introduced in Sub- 
section |^2 which allows the time dependent coupling func- 
tion /i(t, r) to be determined numerically. 

The interaction picture time evolution operator 



t/2B(i,T) 



7-g-i/J<it'H2B(t') 

C/ll(i,T) t/l2(i,T) 



(|bg>,|cl)) 



t/2l(i,T) U22{t,T) I (Cl| 



(bgf 
cl 

(A30) 



corresponding to the two channel Hamiltonian H2B{t) in 
Eq. (0 may be written in terms of the background and closed 
channel Green's functions of the uncoupled motion {W = 0) 
by considering the off diagonal coupling as an inhomgeneous 
term in the Schrodinger equation. The matrix elements rel- 
evant for our purposes are then determined by the following 
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closed set of coupled integral equations: 

C/ii(t,r) = [/bg(i-r) 

+ / dt'G,,g{t-t')W{t')P,e''^'^''^U2i{t',T), 

(A31) 

C/2i(i,r) = -/ dt'e-"^(*')p,W^(i')C/ii(i',T), (A32) 



Here, C4>g(i) = cxp[—iH]^gt/h] is the background chan- 
nel time evolution operator, Ghg{t) — —i9{t)Uhg{t)/h the 
corresponding Green's function of the Schrodinger operator 
ihd/dt — i?bg, and ip{t) — ~ Jf dT{Ei, — hio)/h is the phase 
angle describing the time evolution of the closed channel state 
101,) (cf. Eq. 0). We have furthermore assumed that the de- 
tuning A == oj — JieEy/h of the laser remains constant in 



time. While this assumption is chosen in accordance with 
the experimental sitation of Ref. |8| and simplifies the ex- 
pressions, the general case of varying detuning may be im- 
plemented straightforwardly. 

In terms of the matrix elements Un and U21 the coupling 
function /i(i, r) reads 



h{t,T) = {2TThfe{t- 



(0|l-bg(7bg(i-r)|0) 



dt'^(0|C/bg(i-t')M^(O^.e'^(*')[/2i(i',r)|0) 

(A33) 

To evaluate the integral equations ( IA3H -( lA33> we have used 
the separable ansatz for Vbg and M^|0,y), Eqs. JAll i and JA2l i. 
respectively. We are thus led to the coupled set of integral 
equations 



( unit, r)\_f u^g{t -t)\_ [' , rjT-r^ \d ( Kg{t - t') 

\u2,{t,T))-\ ; X ^'^'Ust'U^ 



Kx{t - 1') 



Ull{t' ,t) 
U2l{t',T) 



(A34) 



for the functions 



in terms of 



un{t,T)=C{x\Uii{t,T)\0), (A35) 

M2i(t,r) =e*'^(*)(<^,|C/2i(t,r)|0), (A36) 



"bgW = CixlUbgim), 



J 



(a/9t)/ibgW = C'(x|Gbg(i)|x>, 



hci{t) 



gi(A+j:7/2)t 

-h{A + ij/2)' 



(A38) 
(A39) 



Inserting the separable ansatz into Eq. ( IA33> the equation for 
(Aj /) jjjg coupling function h{t, t) reads 



d_ 
'dt' 



h{t,T) = Kg{t ~ t) ~ e{t ~ t) I dt'y/l{F)u2lit',T)— e{t ~ t')Ukg{t ~ t') 

I 



(A40) 



where 



V(i) - (27r;i)30(i)(O|FbgC/bg(t)|O>. 



(A41) 



The calculation of the coupling function h{t, t) therefore re- 
duces to evaluating the functions /ibg, Wbg, and /ibg. Note that 
all what is needed of the information included in the full scat- 
tering potential Vbg are matrix elements between the Gaussian 
wave functions constituing the separable ansatz of Eqs. (IA1> 
and JA2> . Hence, the minimal set of parameters introduced 
in Subsection IA ll determines the full dynamics of the low en- 
ergy scattering of a pair of atoms near the photoassociation 
resonance. 

In all many body applications considered in this paper the 
background channel Hamiltonian does not change in time. 
Hence, the functions /ibg, Wbg, and /ibg only depend on a sin- 



gle time variable i — r and may therefore be evaluated us- 
ing the spectral decomposition of the time dependent Green's 
function. 



^AE,\t/h \<ib^){<ib: 



»(pV™)t/?i|^bg^^^bg| 



E' 



(A42) 



Here, \(fi^^) and \(f)^) are the background channel scattering 
and bound states, with eigenenergies p^ /m and Ex, respec- 
tively. Taking only into account the contribution from the scat- 
tering states, i.e. the first term on the right side of Eq. ( IA42> . 
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we obtain 



/ibg(i) - Oit) 



47r?i^abtr 



+ {2nh)HTTmJ^ dp |(O|ybg|0^^)|' exp f-i^ 

(A43) 

The contributions from the discrete part of the spectrum may 
be neglected since all binding energies E\ of •^^Na pairs in 
the (_F = 1, mp = — 1) ground state are much larger than any 
other energy scale introduced through the variation of the laser 
intensity. We therefore introduced a slowly varying envelope 
approximation by consistently neglecting all terms oscillating 
like cxp{i\E^i\t/h} or faster 

The matrix element of Vbg is rewritten as a matrix ele- 
ment of the T-matrix between plane wave momentum states, 
{0\Vhg\(j}l^) = (0|Tbg(pVm)|p). For the low scattering en- 
ergies in a condensate only the partial wave with vanishing 
angular momentum contributes such that this matrix element 
only depends on the modulus of p such that the angular in- 
tegral could be evaluated when deriving Eq. ( IA43> . Finally, 
/ibg {t) was determined by applying the separable parametriza- 
tion introduced in Sect. I A II i.e. Eq. JA7l i for the T-matrix, 
and thus expressing the matrix element in Eq. (IA43> in terms 
of the expectation values in Eq. ( IA6> . 

In an analogous way, we determined the functions Ubg 
and /ibg. Given these functions, uu and U21 were calcu- 
lated iteratively according to the Volterra-type integral equa- 
tion ( IA34> . for a particular time profile /(t) of the laser in- 
tensity. Both functions are defined over the triangular plane 
to < T < t < imax, where <o and t^ax mark the beginning 
and the end of the laser pulse. At the end, h{t, r) is deter- 
mined by use of Eq. ( IA4Qt . 



APPENDIX B: EIGENSTATES OF A NON-HERMITIAN 
OPERATOR 

The two channel Hamiltonian, Eqs. @, Q, describing an un- 
stable molecular state \<j)^) in the closed channel coupled to 
the atoms scattering in the background channel potential Vbg 
is non-Hermitian due to the finite lifetime of \(j>p). Thus the 
corresponding time evolution operator is non-unitary, i.e. the 
norm of any physical eigenstate of the Hamiltonian may de- 
crease in time. 

In this appendix we shall, in Subsect. |BJ recall the basic 
properties of the eigenstates of non-Hermitian operators. In 
Subsect. Ill D 21 we shall discuss in further detail the dressed 
states, i.e. the spectrum of eigenstates of the Hamiltonian in 
Eqs. 0, ^ close to threshold. 

A suitable approximation method to describe the temporal 
behaviour of a number of discrete, partially degenerate states, 
which decay, induced by a weak perturbation, into a range of 
states with a continuous energy spectrum, relates back to the 
work of Weisskopf and Wigner 1 7 1 ] . For not too small and not 
too large evolution times this approximation leads to the usual 



exponential decay of the energy eigenstates, described by a 
decay width contributing as an imaginary part to the energies: 



E, 



:7„ J = 1,2,. 



(Bl) 



The Hamiltonian may then be expressed in terms of its eigen- 
states as 



H = Y, \'l^j)Ej 



(B2) 



where \(j)i) and {(f>i\ are the right and left eigenstates of H to 
the same eigenvalue, respectively, which are defined by the 
relations 



H\q^,) = E,\<j,,), 



\H^ 



h\Ej- 



(B3) 
(B4) 



Note that |0i) is, in general, not transformed into {(t)i\ by a 
Hermitian conjugation. The left and right eigenstates form 
sets of eigenvectors which are dual to each other in the sense 
that they may be orthonormalized by 



i\4>]) = %■ 



(B5) 



Furthermore, the absolute norm of each eigenstate may be 
fixed at a certain time, by requiring 



1. 



(B6) 



While the occupation probability {(f)j{t)\(f)j{t)) exponentially 
decays in time in the Schrodinger picture, the dual norm 
{(f)j{t)\4)j{t)) remains equal to one. Hence, the normaliza- 
tions ( IB5> and ( IB6I 1 give the states a physical meaning and 
provide them with properties which are necessary for basis 
set expansions. For such expansions one introduces the quasi 
projectors 



P,^\<Pi){hl 



(B7) 



with the required properties PiPj = ^ijPj and TrPj = 1. 

In the case that the eigenstates of H with complex energies 
are coupled to each other by a Hermitian operator, the result- 
ing new right energy eigenstates will in general not be orthog- 
onal to each other even if the uncoupled states happened to 
be so. However, a dual basis set of left eigenstates may be 
defined which are orthonormal to the set of right eigenstates. 

In the situation considered in this article, a single unstable 
molecular state is coupled to the spectrum of eigenstates of 
the background Hamiltonian consisting of a discrete subset of 
bound states and a continuous part of scattering states. 



APPENDIX C: MULTICHANNEL FORMALISM OF THE 
MICROSCOPIC MANY BODY THEORY 

Consider a many body ensemble of atoms which may be in 
a defined number of internal states and collide through the cor- 
responding binary scattering channels. We label the different 



internal states by Greek indices. In the case of photoassocia- 
tion in sodium the possible asymptotic internal states will be 
characterized mainly by the electron-nucleus relative angular 
momenta as well as the total atomic angular momenta of the 
collision partners. 

The photoassociation laser is treated classically, and any 



J 
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coupling of an atom to the continuum of photon modes lead- 
ing to spontaneous emission is taken into account, in the 
Wigner-Weisskopf approximation, as a finite width of the rel- 
evant atomic energy levels. The many body Hamiltonian of 
the atoms may be written as 



^W ^ E/ rfxV'Ux)i?^^(x)V.(x) + \Y.J rfxdyV;I(x)V^l(y)F„A,M.(x-y,i)V^4x)V.(x). (CI) 



K.X,fJ,,U 



The operators for creation and annihilation of a single atom 
satisfy the bosonic commutation relations 

[^,(x),V'l(x)]=(5,A'5(3^(x-y), (C2) 

[V'k(x),V'a(x)]-0. (C3) 

The one body Hamiltonian accounts for the kinetic energy of 
an atom, its energy inside an external trap, and the relative 
internal energy relative to some arbitrary zero: 



H 



IB/ 



2m 



+ X4rap(x) 



E, 



int 



(C4) 



The coupling functions Vka,/ji/(x, f) in the interaction term 
in Eq. (ICU contain all information about the incoming and 
outgoing microscopic scattering channels which may be la- 
belled by the pairs of asymptotic internal quantum numbers 
{/UJ/} and {kX}, respectively. For the case of photoassoci- 
ation these are the open background and resonantly coupled 
closed channels as well as the laser induced coupling which 
depend on the laser parameters and the Franck-Condon over- 
laps. The Hamiltonian includes an explicit time dependence 
which enters through the temporally variable laser couplings. 
The theory is formulated in the Schrodinger picture, where 
the density matrix depends on time. All physical quantities 
describing the many body properties are expressed through 



J 



expectation values of normal ordered products of the time in- 
dependent field operators with respect to the state at time t, 
denoted as (•)(. For the purpose of this paper we will con- 
sider, besides the mean field 



*«(x,i) = (V'Jx))t 



(C5) 



the thermal and anomalous density correlation functions 
('0l(y)V'K(x))t and (V'A(y)V'K(x))t, respectively. 

The first characteristic of the theory is a rewriting of the 
infinite hierarchy of non-linear coupled dynamical equations 
resulting for these correlation functions from the Hamiltonian 
dCH in terms of noncommutative cumulants |10, 11|. The 
cumulants to be considered in the following are the mean field 
^k(x, <), Eq. (IC5> . the pair function 

$,a(x, y, t) = {My)M^))t - *.(x, t)*A(y, i), (C6) 
and the density function of the non-condensed fraction: 

r,A(x, y, t) = (7Al(y)^«(x))t - vl/,(x, t)vl/* (y, i), (C7) 

The second characteristic of the approach is a consistent 
truncation of the system of equations llOt Hill . Within the 
leading order approximation which we will apply in this work 
the dynamics is determined solely by the equations 



zn-*„(x,t)=i7iB(x)*,(x,t)+ Y, J dy*l(y,i)F„A,M,(x-y,t)[$^,(x,y,t) + *,,(x,f)*,(x,i)], 
d 

*fi^«'«A (x, y,t) = Y, [^^L'^ (^' y' *)^t-'^ (^' y' *) + ^«^^M^ (x - y, 0*M (x, t)*^ (x, t)] . 

I 



(C8) 
(C9) 



/i,i/ 



where the elements of the Hamiltonian matrix of two interact- 
ing atoms are defined as 



The non-condensed atoms are, to this order, described in terms 



i72B^,(x,y,i) = [i7iB(x)+i7f (y)] +y.A,M.(x-y,t). 

(CIO) 
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of the pair function as: 

r,A(x,y,<) =5] y rfx'$«^(x,x',i)$*^(y,x',t). 



(Cll) 



The total number of atoms at time t. 



Nit) 



E 



J dx [|*,(x, t)|2 + r,,(x, X, t)] (C12) 



is conserved in time if the overall system is closed, i.e. if 
there is no loss due to, e.g., spontaneous photon emission, to 
an atomic level not counted to the system. When consider- 
ing the system of one open background channel coupled to 
a closed channel vibrational state, as described in detail in 
Sect. IllBI we will find that N is, in general, not conserved 
due to the nonzero width 7 of the closed channel state. Dis- 
regarding this loss for a moment, Eqs. ( IC11> and (|^^ reveal 
that the appearance of a non zero pair function, which reflects 
the buildup of pair correlations, means that atoms have been 
lost from the condensate to the non-condensed fraction. These 
correlated atoms become part of a relatively hot, unequili- 
brated cloud around the condensate. In the PA process, where 
7 > 0, they account for the difference between atoms lost 
from the condensate and those transferred into ground state 
molecules through the spontaneous decay process. We finally 
note that inelastic loss phenomena due to e.g. three body re- 
combination fl(?, d(?] are neglected in the leading order ap- 
proach presented here. 



APPENDIX D: TIME EVOLUTION OF THE TOTAL ATOM 
NUMBER 

In the absence of decay the total number of atoms A^tot (t) = 
J dxntot(x; t) in the condensate and the excited states, with 
ntot given by Eq. J38> , is constant in time. This general prop- 
erty of the microscopic dynamical theory was shown for the 
one-channel case in Ref. 1 1 1 ] and generalized for two chan- 
nels in Ref. 14511 . For photoassociation which involves res- 
onant coupling to an unstable closed channel state, however, 
the system is no longer closed and the total number of atoms 



may decrease in time. In this appendix, we sketch the line of 
arguments leading to the evolution equation J39> for the total 
atom number. 

According to Refs. [11,^] the time derivative of the total 
number of non-condensed atoms r(x, y, t) ~ rbg(x, x, t) + 
rci(x, y. t) may be written as 



d_ 
dt 



Nncit) 



d_ 
dt 



= 2Jlm 



/ (ixr(x,x,t) 

dxdy$*(x,y,t) — $(x,y,i) 



(Dl) 



where $ = |bg)$bg + |cl)<I'ci includes the background and 
closed channel pair functions. Using Eqs. J28> and ( I29> for the 
time evolution of the mean field and pair function one finds: 



d_ 
dt 



(7Vc(i) +7v„eW) = mmHMt) - HUt)]mt)), 



(D2) 



where Nc — J dx\'^{x)\^ is the number of condensed atoms 
and |$(i)) denotes in a basis independent way the pair func- 
tion $ (cf. Sect. nil A2t . For obtaining Eq. ( ID2> it is essential 
that there is no direct decay from the background channel, 
i.e. Vj — Vhg, and that a condensate only exists for atoms in 
the internal ground state, as discussed in Sect. IIII ATI 

The term on the right of Eq. ( ID2> is non-zero due to the 
imaginary part of the closed channel diagonal element of the 
Hamiltonian in Eq. ( 13 0> . Hence, inserting Eq. ( I30l l, we obtain 



|(iV,(i)+7V„eW)=-J|(cl,0.|$W)p. 



(D3) 



Here, we have once more used that the hermitian conjugate 
of the time independent right eigenstate \(f>u) is equivalent to 
the left eigenstate ((/)^|. Finally, solving Eq. M9\ one obtains 
an expression for <I>ci(x, y, t) which by insertion into Eq. jD3t 
yields the desired relation (I39> . Thereby, the zero momentum 
plane wave state results from the approximation [11, 45] that 
the mean field ^(x, t) is constant on the scale of the varia- 
tion of the background potential Vhg{r) and the closed chan- 
nel state 4)^{r). 
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